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Nanofluidics has emerged recently in the footsteps of microfluidics, following the quest of scale 
reduction inherent to nanotechnologies. By definition, nanofluidics explores transport phenomena 
of fluids at the nanometer scales. Why is the nanometer scale speciflc ? What fluid properties are 
probed at nanometric scales ? In other words, why 'nanofluidics' deserves its own brand name ? In 
this critical review, we will explore the vast manifold of length scales emerging for the fluid behavior 
at the nanoscales, as well as the associated mechanisms and corresponding applications. We will in 
particular explore the interplay between bulk and interface phenomena. The limit of validity of the 
continuum approaches will be discussed, as well as the numerous surface induced effects occuring 
at these scales, from hydrodynamic slippage to the various electro-kinetic phenomena originating 
from the couplings between hydrodynamics and electrostatics. An enlightening analogy between ion 
transport in nanochannels and transport in doped semi-conductors will be discussed. 



I. NANOFLUIDICS, SURROUNDING THE 
FRAME 

Nanofluidics, the study of fluidic transport at nanome- 
ter scales, has emerged quite recently in the footsteps of 
microfluidics. Pushing further the limits of fluidic down- 
sizing is an attracting stream, in the spirit of scale re- 
duction inherent to all micro- and nano- technologies. 
Various reasons may be seen as to motivate these novel 
developments. First, from the point of view of biotech- 
nological ("lab on a chip") applications, decreasing the 
scales considerably increases the sensitivity of analytic 
techniques, with the ultimate goal of isolating and study- 
ing individual macromolecules [HEI- But also, from the 
point of view of fluidic operations, nanometric scales al- 
low to develop new fluidic functionalities, taking explicit 
beneflt of the predominance of surfaces. Typical exam- 
ples involve preconcentration phenomena [3], the devel- 
opment of nanofluidic transistors [U [S] or the recently 
proposed nano-fluidic diodes [HI [7]. But the analogy to 
micro-electronics is somewhat limited: fluid molecules 
are not electrons and the notion of large scale integration 
for fluidic devices, i.e. nanofluidics as a way of increasing 
the density of fluidic operations on a chip even further, 
is probably not a pertinent goal to reach for fluidic oper- 
ations. 

But from a different perspective, nanofluidics also car- 
ries the hope that new properties will emerge by taking 
benefit of the specific phenomena occurring at the small- 
est scales: new solutions may be obtained from the scales 
where the behavior of matter departs from the common 
expectations. The great efficiency of biological nanopores 
(in terms of permeability or selectivity) is definitely a 
great motivation to foster research in this direction [HUH]- 
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There is indeed a lot of room for improvements at these 
smallest scales: the example of aquaporins (AQP) is in- 
teresting in this context. Aquaporins channels are a key 
component of many biological processes [10 and play the 
role of water filters across biological membranes. These 
channels fuUfill the confiicting tasks of being both ex- 
tremely permeable to water, while extremely selective for 
other species [8]. To give an order of magnitude, the per- 
meability of the water channel is typically 3 orders of 
magnitude larger than what would be expected on the 
basis of the classical fluid framework for the same pore 
size im [T^ . A potential answer is that AQP, although 
a fllter for water, is mainly hydrophobic, i.e. water re- 
pellent ! Of course some hydrophilic polar nodes are dis- 
tributed along the pore, so as to keep water in the mainly 
hydrophobic environment. Quoting the terms of Sui et 
al. in Ref. [B], "i/ie availability of water-binding sites at 
these nodes reduces the energy barrier to water transport 
across this predominantly hydrophobic pathway, while the 
relatively low number of such sites keeps the degree of 
solute-pore interaction to a minimum. In balancing these 
opposing factors the aquaporins are able to transport wa- 
ter selectively while optimizing permeability." Beyond 
this elementary picture, understanding how AQP fuUflUs 
its challenging properties would deflnitely be a source of 
inspiration and open new perspectives for technological 
breakthrough in flltration, desalination, power conver- 
sion, ... However reproducing such a delicate composite 
(patched) architecture in bio- mimetic membranes is a 
great challenge, which requires breakthrough in the con- 
ception of its elementary constituents. But it points out 
that surfaces and their chemical engineering are a key 
actor to optimize fluid properties at nano-scale. 

However, in its roots, nanofluidics is not a new fleld: as 
judiciously recalled by Eijkel and van den Berg in their 
pioneering review on the subject pp, many 'old' fields of 
physics, chemistry and biology already involves the be- 
havior of fluids at the nanoscale. Like Monsieur Jourdain 
in "Le bourgeois gentilhomme" by Moliere, one has done 



2 



'nanofluidics' for more than forty years without know- 
ing it |151j . One may cite for example the domains of 
electro-kinetics (electro-osmosis or -phoresis, ...) with 
applications in chemistry and soil science, membrane sci- 
ence (ultra- filtration, reverse osmosis, fuel cells, ...), col- 
loid chemistry, and of course physiology and the study of 
biological channels 1 . An interesting question is accord- 
ingly wether - on the basis of the novel 'nanofluidic' point 
of view - one may go beyond the traditional knowledge 
in those 'old' fields and obtain unforeseen results, e.g. al- 
lowing for better optimization of existing technologies ? 
Our belief is that the answer to this question is already 
positive and this is one of the key aspects that we shall 
discuss in this review. 

Finally it is also important to note that nanofluidics 
has emerged recently as a scientific field (i.e. naming a 
field as " nanofiuidics" ) also because of the considerable 
progress made over the last two decades in developping 
nano-fabrication technologies, now allowing to fabricate 
specifically designed nanofiuidic devices, as well as the 
great development of new instruments and tools which 
give the possibility to investigate fluid behavior at the 
nanometer scale. One may cite for example: new electri- 
cal detection techniques. Surface Force Apparatus (SFA), 
Atomic Force Microscopy (AFM), nano-Particle Image 
Velocimetry (nano-PIV) coupling PIV to TIRF set-up 
(Total Internal Reflection Fluorescence), as well as the 
considerable progress made in computational techniques, 
like Molecular Dynamics simulations. It is now possible 
to control/design what is occuring at these scales, and 
observe/measure its effects. This is the novelty of the 
flcld, and the reason why nanofluidics now deserves its 
own terminology. 

The paper is organized as follows: In the second sec- 
tion we will replace nanofluidics in the perspective of the 
various length scales at play in fluid dynamics. We shall 
in particular discuss the limits of validity of continuum 
(e.g. hydrodynamic) descriptions. In the third section 
we discuss the dynamics of fluid at interfaces and the 
nanofluidic tools which have been developped recently to 
investigate it. In the fourth section we will explore var- 
ious transport phenomena occuring in diffuse layers. In 
the fifth section we raise the question of thermal noise in 
nanofluidic transport . Finally we will conclude by explor- 
ing some general consideration and expectations about 
nanofluidics, especially in terms of energy conversion and 
desalination. 

As a flnal remark, this review, like any review, is our 
subjective and personal view of the field of nanofluidics 
and the perspectives one may foresee. We organized our 
exploration around the length scales underlying fluid dy- 
namics at the nanometric scales and how nanofluidics 
allows to probe the corresponding mechanisms. Accord- 
ingly, our aim is not to explore exhaustively the - al- 
ready large - litterature of the domain, but merely to 
disentangle the various effects and length scales under- 
lying the behavior of fluids at the nanometer scales. In 
doing so, we certainly hope that this review will raise 



new questions, open new directions and attract people in 
this fascinating domain. 



II. LIMITS OF VALIDITY OF CONTINUUM 
DESCRIPTIONS AND NANOFLUIDIC LENGTH 
SCALES 

The introduction of the terminology "nanofluidics" 
(furthermore to define a specific scientific field) suggests 
that something special should occur for the transport of a 
fluid when it is confined in a channel of nanometric size. 
This leads to an immediate question: why should the 
nanometer length scale have anything specific for fiuidic 
transport ? Nanofluidics "probe" the properties of fluids 
at the nanoscale: so, what does one probe specifically in 
the nanometer range ? 

Actually one may separate two different origins for 
finite-size effects associated with nanometer scales: bulk 
and surface finite-size effects. The former, bulk effects 
are intimately associated with the question of validity 
of the classical continuum framework, in particular the 
Navier-Stokes equations of hydrodynamics: when do such 
descriptions breakdown ? Can one then expect 'exotic' 
fluid effects associated with the molecular nature of the 
fluid ? On the other hand, surface effects play an in- 
creasingly important role as the "surface to volume ra- 
tio" increases {i.e. as the confinement increases). We 
already pointed out the importance of such effects on the 
example of AQP water channels. As we show below, the 
surface effects occur at much larger scales than the 'bulk' 
deviations from continuum expectations. 

The discussion on the length scales at play, to be ex- 
plored in this section, is summarized in Fig. [l] 
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FIG. 1: Various length scales at play in nanofluidics. 



A. Validity of bulk hydrodynamics 

We first start with a short discussion on the validity of 
bulk hydrodynamics. In practice, this raises the question: 
when do the Navier-Stokes (NS) equations break down ? 
These equations were developped in the 19*'' century to 
describe fluid flows at, say, the human scale. What comes 
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as a surprise is their incredible robustness when applied 
to ever smaller scales. 

As a fact, for simple liquids, the continuum frame- 
work of hydrodynamics is apparently valid down 
to the nanometer scale. In other words, there is no 
expected deviations to the bulk NS equations for con- 
finement larger than Inm. This very surprising re- 
sult is actually suggested by a number of experimental 
and molecular simulations studies. On the experimen- 
tal side, one may cite the early work by Chan and Horn 
[13] and later by Georges et al. [U] using Surface Force 
Apparatus, where the prediction of hydrodynamics (the 
Reynolds formula in their case) was verified to be valid 
for confinement larger than typically ten molecular diam- 
eters. More recently, and specifically for water, the works 
by Klein et al. [TS] and E. Riedo et al. [TBI showed that 
water keeps its bulk viscosity down to ^ 1 — 2nm, with 
a drastic change of behavior for stronger confinements, 
where the wettability of the confining surface plays a role 
[16j . A similar behavior was found for other liquids like 
octamethylcyclotetrasiloxane (OMCTS) [iTllTS] . 

This threshold for the applicability of continuum hy- 
drodynamics was also investigated using Molecular Dy- 
namics simulations of confined water, and the same 
value of about 1 nm came out of the simulation results 
[19] [20] [21]. Furthermore, it is interesting to note that 
beyond the validity of continuum equations, the value of 
the viscosity also remains quantitatively equal to its bulk 
value. We show in Fig[2]the results of MD simulations for 
the viscosity 77 of water (SPC-E model) measured in var- 
ious confinements. The value of the viscosity is obtained 
by measuring the shear-stress on the confining plates for 
a given shear-rate (corrected for slippage effects). As 
shown on this figure, the shear viscosity of water keeps 
its bulk value down to confinements of « Inm (typically 
3 water layers). 
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FIG. 2: Viscosity of water versus confinement, as measured 
from MD simulations of water confined in a nano-slit made 
of either two hydrophilic walls, or one hydrophilic and one 
hydrophobic wall (janus confinement), [courtesy of D. Huang] 

Now, in contrast to the viscosity, other transport co- 



efficients may be more strongly affected by confinement. 
This is the case for example of the (self-)diffusion coef- 
ficient which was shown to depend algebraically on the 
confinement width '55', '57. As a consequence, the dif- 
fusion coefficient in confinement strongly departs from 
its Stokes-Einstein prediction, D = kBT/{STT'i]a), with a 
the molecule diameter. The latter is accordingly not a 
correct measure of the viscosity, as sometimes assumed 

mill]. 

Fundamentally the validity of NS equations down to 
typically Inm (for water) is a priori unexpected. Navier- 
Stokes (NS) equations have been developped to account 
for the ffuid dynamics at 'large' scales and they rely on 
an assumption of a local and linear relationship between 
stress and velocity gradients. 

To our knowledge, there is no firmly grounded argu- 
ment for this validity and we propose here a few ten- 
tative leads. NS equations, like any continuum frame- 
work, rely actually on the key assumption of a separa- 
tion of length- and time- scales between the investigated 
length scale and the 'molecular' dynamics. This is the 
"hydrodynamic limit" , in which the continuum frame- 
work should hold. Under this time- and length- scale 
separation, the microscopic dynamics associated with a 
huge number of degrees of freedom Af ~ 10^^ reduces 
to equations with just a few degrees of freedom (veloc- 
ity field, pressure, density, temperature, etc.), while all 
the complexity is hidden in just a few phenomenologi- 
cal coefficients. Achieving this averaging out of the "fast 
variables" is a huge challenge which can be achieved sys- 
tematically only in a few limiting cases, see e.g. |25| 
for an explicit example. In general, the elimination of 
fast variables is summarized in the Green-Kubo relation- 
ship for the transport coefficients: for example the shear 
viscosity r] can be expressed as the integral of the stress- 
stress correlation function according to 

1 f°° 

V = Yk^x X (f^xyW^^a(0)>oqudi (1) 

with axy a non diagonal component of the stress tensor 
[2S] and (•)equ denotes an equilibrium average. 

The time-scale separation underlying the validity of 
Eq. ^ requires the microscopic time scales. To-, charac- 
terizing the stress-stress correlation function, to be much 
smaller than a hydrodynamic time scale. This time scale 
is e.g. the relaxation time of momentum, which for a 
given wave vector q has the expression Tg = (vq^)^^, 
with V — "q/ p the kinematic viscosity and p the mass 
density. An implicit condition of validity is accordingly 
vq^T„ < 1. This fixes the limit for time-scale separation 
at confinements w larger than a viscous length scale: 

u; > 4 = • (2) 

Putting numbers for water, a typical correlation time for 
the stress-stress correlation function is in the picosecond 
range To- ~ 10~^^s, while ly = 10^^ m^.s^^. Altogether 
this gives for water 

ia ~ Inm (3) 
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A nanometric characteristic length scale thus emerges 
naturally as a lower bound for the validity of the notion 
of viscosity, and thus for the apphcation of standard NS 
hydrodynamics. This is indeed the experimental value 
below which strong departure from the continuum frame- 
work is observed. 

Of course, it would be interesting to explore further 
the physical contents of this condition: what occurs for 
smaller confinements ? Also, can one tune Ic in specific 
systems ? etc. 

In the extreme limit of single file transport, the trans- 
port behavior is indeed expected to strongly deviate from 
the bulk expectations, with the predicted occurence of 
anomalous diffusion, and a fast stochastic transport, to 
cite some examples [28l |29l [30]. Apart from biological 
systems -where molecular pores are the rule fabricat- 
ing channels with molecular size is however still out of 
reach with state-of-the-art nanofabrication technologies. 

As a side remark, it would be interesting to extend 
this discussion to other transport phenomena, like for 
example heat transfer at nanoscales. Heat transfer from 
nanoparticles and nanostructures to a fluid is an actively 
growing field, e.g. in the context of cooling enhance- 
ment by structured surfaces or local heating of fluids via 
nanoparticles |27j . Similar arguments as above may be 
proposed, and lead to an equivalent sub-nanometer heat 
length scale (for water at room temperature). This sug- 
gests that the continuum phenomenological picture for 
heat transport in fluids, involving Fourier law, is expected 
to hold down to nanoscales, in agreement with observa- 
tions [27] . However a full investigation of the limit of ap- 
plicability of the corresponding continuum phenomeno- 
logical laws remain to be developped in general. 

A conclusion of the above discussion is anyway that for 
most nanofluidic applications involving water in supra- 
nanometric confinements, bulk NS equations can be 
safely used to account for the fluid transport. 



B. A broad spectrum of length scales 

Aside from introduced above, there is a number of 
length scales which enter nanofluidic transport. These 
length scales are all related to surface effects in one way or 
another. We review in this section a panel of these length 
scales, climbing up from the smallest - molecular - scales 
to the largest, but still nanometric, scales. This panel is 
not exhaustive, and we merely insist on the main length 
scales which appear generically in nanofluidic problems. 
Our aim is to show that there are a number of pertinent 
scales which lie indeed in the nanometric range. This 
implies that speciflc nanofluidic phenomena will show up 
when the conflnement compares with these values. 



1. Molecular length scales 

At the smallest scales, the granularity of the fluid 
and its components (solvent, ions, dissolved species, ...) 
should play a role. This is defined by the molecular scale, 
associated with the diameter of the molecules a, typically 
in the angstrom scale, SA for water. As we discussed 
above, this is the scale where the validity of hydrody- 
namics breaks down, typically for confinements of 3-5 
molecular diameter. Various phenomena then come into 
play. First, the structuration and ordering of the fluid at 
the conflning walls plays a role. Such ordering was indeed 
shown experimentally to induce oscillatory dissipation in 
liquid films with a width of several molecular diameter 

[ni[i8]. 

Similarly non-local rheological properties have been 
predicted in strongly confined situations [3T], associated 
with a length characterizing non-locality typically in a 
few diameters range (indeed in line with the above limit 
of NS hydrodynamics) . This result echoes recent findings 
of non-local fiow curves in confined soft glasses [S2]. An 
interesting consequence of non-locality, as pointed out in 
[32] , is that the specific nature of surfaces does influence 
the global flow behavior in confinement, in line with re- 
cent observations in AFM measurements by Riedo et al. 
for strongly confined water (less than Inm) films at a 
hydrophilic and hydrophobic surfaces [TfT . 

But even more drastic are the deviations occuring in 
the extreme limit where only one molecule can enter 
the confining pore. Very strong correlations and col- 
lective motion builds up in the liquid dynamics, lead- 
ing to the so-called "single-file" transport [33J[51]. The 
associated transport differs strongly from the bulk hy- 
drodynamic predictions and various specific phenomena 
show up, like non-fickian transport [34 , fast water trans- 
port and ion transfer, stochastic fiow with the occurence 
of "bursts", etc., see e.g. [33,- Simulations of water 
transport through single-file carbon nanotubes (where 
only one water molecular can pass through the channel) 
showed indeed flow rate much larger than those predicted 
from the Poiseuille law of hydrodynamics [28] . 

Note however that this behavior is limited to 
nanochannels where molecules cannot cross each other, 
i.e. with diameters not exceeding the water molecular 
length scale, like in biological channels. Up to now, ar- 
tiflcially produced nanochannels do not reach this limit. 
This point is interesting to discuss in the context of re- 
cent experimental results obtained for flow through car- 
bon nanotube membrane [Hi 1351 I36j . A recent exper- 
imental flrst was achieved in the group of O. Bakajin, 
who investigated the flow through a membrane consti- 
tuted of carbon nanotube with a size in the range 1.3-2 
nm [9^, following the work on larger nanotubes (~ 7nm) 
by Majumder et al. [33]. These works demonstrated 
massively enhanced flow permeability, as compared to 
bulk predictions (in term of slippage, to be discussed be- 
low, this would correspond to slip lengths up the micron, 
far above any expected result for these surfaces). More 
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recent experiments performed with membranes made of 
wider carbon nanotubes suggest a large, but still smaller, 
enhancement of the permeability |36j . Therefore, while 
it would be tempting to discuss these results in the con- 
text of single-file transport and the predictions by Hum- 
mer et at on fast single-file transport in sub-nanometer 
carbon nanotube [33^, this point of view is difficult to 
justify: as discussed above, fluid flow in nanotube pipes 
with a supra-nanometric diameter can be safely described 
by continuum approaches (Stokes equation) and the dy- 
namics is definitely not single file. Molecular Dynamics 
simulations of flow through similar carbon nanopipes in- 
deed suggest large flow enhancement [211 EZ], but still 
far below the experimental results. There is therefore a 
strong need to perform flow experiments through a single 
object, i.e. through a single carbon nanotube. This is 
still an experimental challenge. 

Altogether transport at the molecular scale offers a 
rich panel of transport behavior. This is definitely the 
scale where, as we quoted in the introduction, new so- 
lutions and technological breakthrough could emerge as 
the behavior of matter departs from the common expec- 
tations. This suggests that the molecular scale would 
be the ultimate scale to reach for nanofluidics. But at 
present, producing channels of sub-nanometric size re- 
mains a technological challenge. 



2. Electrostatic length scales and their dynamtc influence 

The scales associated with the interaction of the flow- 
ing matter with its environement also enter the game. 
And of course, electrostatic plays a key role in this con- 
text. We leave aside interactions deriving indirectly 
from electrostatics to focus on the interactions between 
charged species. These obviously concern ions dissolved 
in the fluid: ion-ion interaction, and interactions of ions 
with charged surfaces. This also involves the dipolar in- 
teraction among water molecules, as well as H-bonding. 

Electrostatics involves a rich panel of length scales. 
And more importantly, due to the long range nature of 
electrostatics, these length scales are intimately coupled: 
effects occuring at the smallest electrostatic lengths do 
climb the scales to span the whole range. 

a. Bjerrum length (and some derivatives) 

We start with the Bjerrum length, £b, which is defined 
as the distance at which the electrostatic interaction Vci 
between two charged species becomes of the order of the 
thermal energy, fc^T: Vci(^b) ~ ksT. For two ions, with 
valence z, embeded in a dielectric medium with dielectric 
constant e, this takes the form 



(4) 



The Bjerrum length therefore compares with the molecu- 
lar range in this case. Note however that for multi-valent 
ions or for organic solvent with lower dielectric constants, 
£ B niay be much larger and thus separate from the molec- 
ular range. 

Other length scales could be introduced along the same 
lines by comparing thermal energy to e.g. charge-dipole, 
dipole-dipole, etc. interactions, taking into account the 
multi-pole interaction under consideration. This leads a 
priori to length scales smaller than the Bjerrum length 
introduced above. For the charge-dipole interaction, this 

1/2 

gives for example £4 — (^sf) , with p the dipole 
strength: for water parameters, £d is in the Angstrom 
range. It does correspond to the typical thickness of a 
hydration layer around an ion |38j . 

By definition, the Bjerrum length (and its derivatives 
like £d above) is the scale below which direct electro- 
static interactions dominate over thermal effects. Its 
consequences on nanofluidic transport is thus expected 
at molecular scale. For example, for confinement below 
£b one expects a large free-energy cost to undress an ion 
from its hydration layer and make it enter a molecular 
pore. This has therefore immediate consequences on the 
filtering process of charged species, as in biological ion 
channels. 

But as pointed out above, a strong interplay with up- 
per scales is expected due to the long range nature of 
electrostatic interactions. For example, as we will dis- 



cuss in section IV ion-specific effects at interfaces, the 
origin of which occurs at a scale £b, have strong effects 
on the electro-kinetics at nanometric, micrometric and 
even larger scales. 

h. Debye length scale 

A central concept of electrostatics is the notion of elec- 
trostatic diffuse layer (EDL) . At a charged interface, the 
EDL is the region where the surface charge is balanced 
by the cloud of counterions. In this region the ion con- 
centration profiles depart from their bulk values due to 
the interaction of ions with the surface charge. This is 
therefore the region where local electro-neutrality is not 
obeyed. 

The Debye length emerges naturally from the Poisson- 
Boltzmann theory and characterizes the electrostatic 
screening in the bulk electrolyte [39]. It is defined in 
terms of the salt concentration according to the ex- 
pression: 



.-1/2 



(5) 



(e the elementary charge). For bulk water at ambiant 
temperature and a valency z = I, this gives £b = 0.7nm. 



Note that in this expression, as in the following, ps de- 
note the salt concentration in the reservoirs, thereby 
fixed by the chemical potential of the ions. The value 
of the Debye length depends on the salt concentration 
and ranges typically between tens of nanometer (30 nm 
for ps = 10~^M) down to Angstroms (3 A for ps = IM). 
In physical terms, the electrostatic free-energy of an ion 
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dressed by its spherical cloud of counterions with size A d 
is of the order fc^T. 

The Debye length characterizes the width of the EDL. 
Interestingly, it is independent of the surface charge and 
only depend on the bulk ion concentration. However, in 
the limiting case of a salt free solution (with Ad oo), 
the width of the diffuse layer is merely fixed by the so- 
called Gouy-Chapman length, to be defined in the next 
paragraph. 

The Debye length plays a central role in nanofluidics 
for various reasons. First, as quoted above, this is the 
region in the fluid close to charged surfaces where local 
charge electroneutrality is broken. Under an applied elec- 
tric field, this is thus the region where volume electric 
forces, /e = PeEe, will apply, with the charge den- 
sity and Ef. an applied electric field. On the other hand 
these forces vanish in the bulk of the material due to 
a vanishing charge density pe ^ far from the surfaces. 
Therefore tuning the fiuid properties inside the nanomet- 
ric EDL is expected to affect the whole response of the 
system under an electric field. Accordingly, tuning the 
dynamics of the fluid in the nanometric EDL, or mod- 
ifying its structure, will have a macroscopic impact on 
the fluid dynamics at scales much larger than the Debye 
length. We will come back extensively on this question 
below, in section [IVB[ 

Second, a specific behavior is expected when the De- 
bye layers overlap in a nanopore, which occurs when its 
size is of the order of twice the Debye length. This phe- 
nomenon has a strong effect on fluidic transport and has 
been the object of intense research recently, as reviewed 
recently by Schoch et al. 
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pO;, see Sec. 

cite in particular the phenomena of ion enrichement and 
exclusion [5]. While these phenomena were known for a 
long while in membrane technology, they found new ap- 
plications in the field of nanofiuidics. It is at the root of 
a number of novel fiuidic phenomena, such as permselec- 
tivity PI , nanofiuidic diodes [B] or surface dominated 
ion transport [41J. 

A key practical aspect is that the condition of Debye 
layer overlap is much less stringent to fuUfiU with actual 
nanofabrication techniques. These are now able to pro- 
duce individual nanopores with at least one dimension in 
the range of ten nanometers, i.e. potentially suited for 
Debye layer overlap. 



centration. Physically, the Gouy-Chapman length can be 
defined as the distance from the wall where the electro- 
static interaction of a single ion with the wall becomes 
of the order of the thermal energy. For typical surface 
charges, say Se ^ 50 mC/m^ 0.3e/nm^) to fix ideas, 
then Hqc ^ Inm. The Gouy-Chapmann length plays a 
role for solutions with very small salt concentration. This 
has an infiuence on the conductance in nanochannels as 
we discuss in Sec. IIVCI 

More interesting is a length which can be defined on the 
basis of the comparison between the bulk to the surface 
electric conductance (relating electric current to an ap- 
plied electric field). This introduces what can be termed 
as a "Dukhin length" , by analogy to the Dukhin number 
usually introduced for colloids |152j . Indeed the conduc- 
tion probes the number of free carriers (ions), so that in 
a channel of width h and surface charge S, the equiva- 
lent bulk concentration of counterions is 2S/ft,. One may 
define a Dukhin number Du — \Yi\/ (hps), where pa is 
the concentration of the salt reservoir |32]- A 'Dukhin 
length' can then be defined as 



M 
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To put numbers, for a surface with a surface charge den- 
sity eS ^ 50mC/m^ (« 0.3e/nm^), iou is typically 0.5 
nm for p^ = IM, while = 5 pva for p^ = 10~''M ! 
The Dukhin length characterizes the channel scale below 
which surface conduction dominates over the bulk one. 
In a different context of charge discontinuities at surfaces, 
it has also been interpreted in terms of an electrokinetic 
"healing" length 

This length scale can actually be rewritten in terms 
of the Debye and Gouy-Chapman lengths as tou ^ 
Xj^/lcc- Note also that in the limit where the Debye 
length is large compared to the Gouy-Chapmann length, 
£gc < ^D, the non-linear Poisson-Boltzmann expres- 
sion for the electrostatic potential [39 allows to rewrite 
the above length as £du ~ Ad exp[e|ys|/2fcsr], with Vg 
the surface potential. The corresponding form for the 
Dukhin number is more standard in colloid science [42'. 

This length plays an important role for the conduc- 
tance in nanochannels, cf Sec jlV C| where surface effects 
are shown to strongly affect conductance IIH |H] . 



c. Length scales associated with surface charge: 
Gouy-Chapman and Dukhin 

Other electrostatic length scales may be constructed 
on the basis of surface electrostatic properties. We al- 
ready mentioned the Gouy-Chapman length scale, which 
shows up for the behavior of a salt solution very close to 
a charged surface. For a surface charge density Se, this 
has an expression £gc — 1/(27''S^b). In contrast to the 
Debye length, the Gouy-Chapman length depends explic- 
itly on the surface density S of charges of the confining 
surface (S in units of m~^), but not on the bulk ion con- 



3. Slip lengths and surface friction 

Up to now we merely considered length scales char- 
acterizing the structure of the fluid and its components 
(ions, ...). However the dynamics of fluids at interfaces 
introduce various length scales. This is in particular the 
case of the so-called slip length, b, characterizing the hy- 
drodynamic boundary condition of a fluid at its confining 
interfaces. The latter is defined according to the Navier 
boundary condition (BC) as [21 

bd„vt=vt (7) 
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with n, t points to the normal and tangential directions 
of the surface, vt is the tangential velocity field and b 
the slip length. The slip length characterizes the friction 
of the fluid at the interface and large slip lengths are 
associated with low liquid-solid friction. 

This point has been amply explored experimentally 
and theoretically: a key result which emerges from this 
measurements is that the slip length of water at solid sur- 
faces depends crucially on the wettability of the surface 
[5T1 US]. We will come back more exhaustively on this 
point in Sec. At this level, one should keep in mind 
that slip lengths in the range of a few ten nanometers 
are typically measured on hydrophobic surfaces, while b 
is sub-nanometric on hydrophilic surfaces. Note however 
that very large slip length, in the micron range, may 
be obtained on nano- and micro- structured interfaces 

gSlIlT]. 

This offers the possibility to modify the nanofluidics in 
pores using chemical engineering of the surfaces. As the 
pore size compares with b a considerable enhancement 
of fluid transport is accordingly expected. Furthermore 
one may remark that these values of b also compares 
with typical Debye lengths. Therefore shppage effects 
are expected to affect ion transport at charged surfaces. 
We will discuss these aspects more exhaustively below. 

A final remark is that other surface related lengths 
may be constructed. This is the case in particular of the 
Kapitza length, which characterizes the boundary con- 
dition for thermal transport across an interface, see e.g. 
[2lJ. The Kaptiza length behaves quite similarly to the 
slip length and strongly depends on the wettability of the 
interface 1481 . 



4. Other length scales 

The review of length scales above is not exhaustive and 
many other length scales may enter the game, depending 
on the problem, geometry and system considered. For 
example, in the context of transport of macromolecules 
(colloids, polymers, DNA, RNA, proteins, ...), the typi- 
cal size of the particle - diameter or radius or gyration - 
plays of course a central role. More precisely a key quan- 
tity is the free energy associated with the confinement 
of the macromolecule, J^c which fixes the partitioning of 
the molecule in the nanopore w.r.t. the bulk. The prob- 
abihty of passage of a macromolecule through a pore is 
expected to scale like exp[— ^c/fcs?"], hence fixing its per- 
meability [49' . For example for a polymer chain confined 
in a pore with size D, the entropic cost of confinement 
takes the form J'c ~ kBT{Rg/ Df/" , with Rg = a x N" 
the radius of gyration of the polymer (a the monomer 
size and N the polymer length). For a polymer in good 
solution, the Flory exponent is ly = 5/3. Pores the size of 
which is below the radius of gyration of the molecule thus 
act as molecular sieve. This has been explored in various 
experiments of polymer translocation in nanopore follow- 
ing the pioneering work of Bezrukov et al. |50l [Flj and 



Kasianowicz[S5], as well as in the context of DNA sepa- 
ration 1531 . 



C. Some practical conclusions 

The discussion above is summarized in Figjl] where we 
organized the various scales at play along the scale axis. 
This points to several conclusions. 

1 . A number of specific phenomena occur in the nano- 
metric range, which indeed justifies the specificity 
of the dynamics of fiuids in the nanometer range, 
i.e. nanofiuidics. 

2. An important conclusion is : For water under 
normal conditions, the Navier-Stokes equation re- 
mains valid in nano-channels down to typically 1-2 
nanometers. This is good and bad news. Good 
because one may safely use NS equations for most 
nanofluidic phenomena. This robustness of contin- 
uum hydrodynamics is remarkable. 

But this is in some sense bad news since it means 
that reaching specific effects associated with the 
granularity of the fluid requires molecular conflne- 
ment, i.e. confinement below the nanometer. Pro- 
ducing nanochannels with such sizes is technolog- 
ically diflicult to achieve. On the other hand, 
this probably explains why most biological channels 
have a molecular entanglement: to reach specificity 
and avoid the "universality" associated with con- 
tinuum equations, one should reach the fiuid molec- 
ular scale. 

3. This suggests that the molecular pore is the ulti- 
mate scale to reach for nanofiuidics, where, as we 
quoted in the introduction, new solutions and tech- 
nological breakthrough could emerge as the behav- 
ior of matter departs from the common expecta- 
tions. But at present, there is however still a long 
way before this scale is technologically accessible. 

4. While it remains difficult to tune the bulk behavior 
of the confined fiuid, there are much more room and 
possibilities to benefit from surface effects. This is 
clearly apparent in Fig. [T| where surface effects 
enter the nanofluidic game at much larger length 
scales. Surface effects enter in particular via the 
Debye length, Dukhin length, slip length, all typi- 
cally in ten nanometers range (and even more for 
the Dukhin length). This suggests that speciflc ef- 
fects will show up when one of these lengths will 
compare with the pore width. Furthermore, par- 
ticular effects should also occur when two of these 
lengths become comparable, independently of the 
confinement. This will be explored in Sec. |IV[ 

Altogether nanofiuidics is a incredible playground 
to play with surface effects ! 
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III. DYNAMICS AT SURFACES AND THE 
TOOLBOX OF NANOFLUIDICS 

We have discussed above that the bulk laws of hy- 
drodynamics are valid down to very small scales, typi- 
cally ^ Inm, so that bulk Navier-Stokes equations can 
be used for most nanofluidic flows. However as the size 
of the nanochannel decreases, the dynamics at its sur- 
face should play an increasingly important role. Navier- 
Stokes equations require boundary conditions (BC) for 
the hydrodynamic flow at the device's surface, and a spe- 
cific knowledge of these BC is a pre-requisite to appre- 
hend fiow at the nanoscale. 

As we quoted above, the BC at a solid surface intro- 
duces a new length, the so-called Navier length or slip 
length, which relates the tangential velocity vt to the 
shear rate at the wall [54]: 



Vt 



dz 



(8) 



Obviously the control of the slip length is of major im- 
portance for flows in confined geometries, and it can have 
dramatic consequences on the pressure drop, electrical, 
and diffusive transport through nanochannels. For in- 
stance the relative increase in hydraulic conductance of 
a cylinder of radius r due to wall slippage is 1 -I- 66/r; 
it can be much larger for the electrical conductance or 
electro-osmosis, as we will discuss in the sections below. 

As a consequence a sustained interest has been devoted 
in the last ten years for the investigation of the BC and its 
dependence on interfacial properties, such as the surface 
topography and the liquid-solid interactions. Reviews on 
this subject are [211 [55| l56]. 



V(Z) 



FIG. 3: Hydrodynamic slip and definition of the slip length h: 
the slip length is the distance in the solid at which the linear 
extrapolation of the velocity profile vanishes. 



A. Theoretical expectations for slippage 

The theoretical understanding of slippage has been the 
object of intense research over the last decade, see Ref. 
[2T] for a review. The BC has been studied theoretically 
by Molecular Dynamics simulations, as well as in the con- 
text of hnear response theory ^451 IFfl [551 [Ml [5U1 [CT] 155] . 



It is now well understood that the BC on atomically 
smooth surfaces depends essentially on the structure of 
the hquid at the interface, itself determined by its inter- 
actions, its commensurability with the solid phase, and 
its global density [21j. 

In order to connect the slip length to the interfacial 
properties, it is useful to interprete the slip length in 
terms of liquid-solid friction at the interface. The friction 
force at the liquid-solid interface is linear in slip velocity 

Vt- 



-A\vt 



(9) 



with A the liquid-solid friction coefficient and A the lat- 
eral area. By definition the slip length is related to the 
friction coefficient according to h — "q/X, with 77 the bulk 
viscosity. As a phenomenological friction coefficient, A 
can be expressed in terms of equilibrium properties in 
the form of a Green-Kubo (GK) relationship [2T1 l59l [62] : 



A 



1 



AksT 



(10) 



where Ff is the total (microscopic) lateral force acting 
on the surface. In practice, this expression is difficult to 
estimate. However some useful information on the slip 
length may be extracted from it. 

A rough estimate of slippage effects can be proposed 
[63] . One may indeed approximate the GK expres- 
sion as A ~ AksT (-^/)cqu ^ with t a typical relax- 
ation time and (i^^)oqu the mean-squared lateral surface 
force. We then write t ~ a^/D where cr is a micro- 
scopic characteristic length scale and D the fluid dif- 
fusion coefficient, while the rms force is estimated as 
(i^j)equ ^ C±p(T{e/a)'^ , with e a typical fluid-solid molec- 
ular energy, p the fluid density and C± a geometric factor 
that accounts for roughness effects at the atomic level 
(large C± corresponding to a larger atomic roughness) 
[45] l63] . Altogether this provides a microscopic estimate 
for the slip length as 



kgTr^D 



(11) 



A more systematic derivation leads to a very similar 
result [2TI [62], with Eq.([Tl|) multiplied by the inverse 
of the structure factor of the fluid, Sw{q\\), computed 
at a characteristic wave- vector of the solid surface q||: 
b oc ^^^(gii)"^. This term measures a kind of commen- 
surability of the fluid with the underlying solid struc- 
ture, in full analogy with solid-on-solid friction. Note 
also that the slip length depends on the product rj x D, 
so that according to the Stokes-Einstein relationship, the 
slip length is not expected to depend on the bulk fluid 
viscosity (except for specific situations where these two 
quantities can be decorrelated) . The scalings proposed 
by the above simple estimate are in good agreement with 
molecular dynamics results [2TJ [63] . 

According to Eq. (Ill significant slippage should oc- 
cur on very smooth surfaces, in case of a low density 
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FIG. 4: Slip length b of water as a function of the contact 
angle on various smooth surfaces. Top : Results from Molec- 
ular Dynamics simulations by Huang et al. for the slip length 
of water on various surfaces, from Ref. [45] • The dashed 
line is a theoretical fit scaling like 6 cx (1 -I- cos6)~^. Bottom 
Compilation of various experimental results. 



at wall - which requires both a moderate pressure and 
unfavorable liquid- solid interactions -, or a low value of 
Sw{q\\), which characterizes liquid/solid commensurabil- 
ity. It also suggests that low energy surfaces, i.e. with 
small liquid-solid interaction e, should exhibit large shp- 
page. Accordingly, hydrophobic surfaces should exhibit 
larger slip length than hydrophilic ones. 

We plot in Fig. |4] (top) results for the slip length 
obtained by Molecular Dynamics of a water model over 
a broad variety of surfaces jlSj [63] . A 'quasi-universal' 
dependence of the slip length on the contact angle is mea- 
sured, with - as expected from the above arguments - an 
increase of the slip length for more hydrophobic surfaces. 



B. Experimental results 

On the experimental side a race has engaged to- 
ward the quantitative characterization of BCs within the 
nanometer resolution, using the most recent developm- 
ments in optics and scanning probe techniques |551 156j . 
This experimental challenge has generated considerable 
progress in the experimental tools for nanofluidic mea- 
surements, and the values of the slip lengths reported 
had a tendancy to decrease as the resolution and robust- 
ness of the techniques has increased. An intrinsic diffi- 
culty stems from the fact that the BC is a continuous 



medium concept involving an hydrodynamic velocity av- 
eraged on many molecules in the liquid and extrapolated 
to the wall: it does not simply reduce to the velocity close 
to the surface of molecules or tracer particules, which 
are furthermore submitted to brownian motion and re- 
versible adsorption. Also, the surface has to be smooth 
and homogeneous on the scale of the probed area in order 
to characterize an intrinsic BC, and not an effective slip 
length. Thus, the problem of the BC has been a discrim- 
inating test of the capacity of instruments to perform 
quantitative mechanical measurements at the nanoscale. 
We propose here a brief overview of the current state-of- 
the-art. This offers the opportunity to quote the various 
tools which have been developped recently to study flow 
and fluid dynamics at the nanofluidic scales. 

The Surface Force Apparatus provided the first sub- 
nanometer resolution measurement of BCs, based on the 
viscous force acting between two crossed cylinders 
whose distance D is measured within the Angstrom reso- 
lution by the so-called FECO interferometric fringes |64j : 



^ 67^r^R2I) 
D 



/*PA)=3^[(l + ^)ln(l + ^)-l)] 



65' 



D' 



(12) 
(13) 



with J] the liquid viscosity and i? the cylinder radius |65j . 
In case of a no-slip BC, reduces to the Reynolds force 
between a sphere and a plane. Chan and Horn [13 in- 
vestigated the flow of various organic liquids conflned 
between atomically flat mica surfaces and showed that 
they obeyed the no-slip BC, with the no-slip plane lo- 
cated about one molecular size inside the liquid phase. 
Their results were confirmed by other experiments run- 
ning the SFA in a dynamic mode, and extended to water 
and to various (wetting) surfaces |14l 166] . Further inves- 
tigations focused on non-wetting surfaces and reported 
the existence of very large slip lengths (~ 1/xm) asso- 
ciated to a complex behaviour, with slippage appearing 
only under confinement and above a critical shear stress 
[67]. However it was shown that elasto- hydrodynamics 
effects due to non perfectly rigid surfaces (such as glued 
mica thin films) retrieve this typical behaviour |68,. Ex- 
periments on rigid surfaces later reported the existence 
of substantial slippage of water on smooth homogeneous 
non- wetting surfaces, with a slip length increasing with 
the contact angle and reaching a range of 20nm at contact 
angles of 110" [Ml EO], see Fig. [i] 

The AFM with a colloidal probe operates along a sim- 
ilar principle but with a smaller probed area than the 
SFA. In pioneering investigations Craig et al. found a 
shear-rate dependent BC of sucrose solutions on par- 
tially wetting SAMs, with slip lengths reaching 20nm at 
high velocities ^71J, while Bonacurso et al. found a con- 
stant (no shear dependent) slip length of 8-9nm for water 
on wetting mica and glass surfaces [72]. These dissimi- 
lar findings also did not agree with available theoretical 
results. However Ducker et al, using an independant 
optical measurement of the probe-distance in the AFM, 
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found a regular no-slip BC on a similar system as the one 
investigated by Craig et al. [73]. It was recently shown 
by Craig et al. that the lateral tilt of the cantilever was 
a significant issue in interpretating hydrodynamic forces 
|74j . Vinogradova et al. also evidenced the important 
effect of the dissipation due to the cantilever [T^ • Using 
a "snow-man" probe made of two spheres glued on the 
top of each other, they showed a perfect agreement with 
the no-slip BC on hydrophilic surfaces, confirming SFA 
experiments, and a moderate slip length of lOnm of water 
onto rough hydrophobic surfaces. Finally dynamic mea- 
surements using an oscillatory drive were implemented in 
the AFM by Maali et al. so as to measure the dissipation 
coefficient with high accuracy [76j. Using this technique 
they found a no-slip BC of water onto mica surface, and 
a 8nm slip length on atomically smooth HOPG graphite 
with contact angle 70° [77] . 

An alternative way to measure interfacial hydrody- 
namic dissipation is through the dissipation-ffuctuation 
theorem, for instance by measuring the modification of 
the brownian motion of a probe at the vicinity of a sur- 
face. In Joly et al. experiment [22] the average transverse 
diffusion coefficient of colloidal particles confined in a liq- 
uid slab was measured by Fluorescence Correlation Spec- 
troscopy as a function of the film thickness, and related 
to the boundary condition on the slab walls. A resolution 
better than 5nm is obtained with particles of 200nm nom- 
inal diameter, by averaging the data over 100000 events 
crossing the measurement volume. The originality of the 
method is to give access to the boundary condition with- 
out any flow, i.e. really at thermodynamic equilibrium. 
It is thus directly comparable to theoretical predictions. 
A collateral result is the direct probe of a property of 
great interest in nanofluidics, i.e. the motion of a solute 
close to a wall. 

Besides the dissipative methods, great progress was 
achieved in the development of optical methods giving 
direct access to the flow fleld and its extrapolation to the 
solid wall with a sub- micrometer resolution. 

The Micro-Particle Image Velocity (/i-PIV) technique 
has been widely used in microfluidics. First attempts 
to determine the BC yielded slip lengths of micrometer 
amplitude [7S] , difficult to conciliate with theoretical 
expectations, although the variation of slippage with the 
liquid- solid interaction were consistent. Gaining better 
resolution raised many issues. Besides colloidal forces 
acting on particles close to a wall, which have to be taken 
into account, such as DLVO forces or electro-osmosis in- 
duced by the streaming potential, an important source of 
uncertainty is the determination of the wall position: lo- 
cating the latter from the signal of moving particles lead 
to a systematic bias because of the depletion induced by 
the colloidal lift [79j . By using the signal of non-moving 
particles adsorbed to the wall Joseph et al [SD] increased 
the resolution on the slip length to lOOnm in a configura- 
tion similar to [78], and found no slippage at this (lOOnm) 
scale whatever the wettability of the surface. A recent 
investigation using double-focus FCS confirmed the exis- 



tence of slippage of a few ten nanometers on hydrophobic 
surfaces [81]. A major improvement toward nanometric 
resolution is the Total Internal Reflection Fluorescence 
method which uses an evanescent wave at the interface 
to be probed, and allow to measure the distance of a 
tracer particle to the wall through the exponential decay 
of the signal amplitude. Bouzigues et al. found a wa- 
ter slip length of 21nm (±12nm) on OTS coated glass 
capillaries [101] . while Lasne et al. measured a water 
slip length of 45nm (±15nm) on a similar surface |82| . 
These results are consistent and the different values of 
slip length may well be ascribed to slight variations in 
surface preparation. Note that in such measurements, 
the size of the probe itself is a major issue in the race to- 
ward nanometric resolution: Brownian motion increases 
as the particle size decreases, not only adding noise to 
the measurement, but allowing particles to leave the fo- 
cal plane between two successive images which decreases 
the locality of the measurement. Some specific statistical 
treatment has to be performed in order to avoid any bias 
in the estimated velocity [82 ] l83 l ITOT] . 

A summary of the results obtained with water on var- 
ious smooth surfaces is gathered in Fig|4] plotting the 
slip length as a function of the contact angle. The var- 
ious methods show the same trend: water does not slip 
on hydrophilic surfaces, and develops significant slippage 
only on strongly hydrophobic surfaces. The highest slip 
lengths reach the range of 20-30nm: accordingly hydro- 
dynamic slippage is non relevant in microfiuidic devices, 
but of major importance in nanochannels. On this ba- 
sis one expects that water flow in highly confined hy- 
drophobic pores, such as mesoporous media or biolog- 
ical channels (c./.the Aquaporin example in the intro- 
duction), should be much less dissipative than ordinary 
Poiseuille flow. The results obtained with the various 
methods show a very good qualitative agreement with 
theory and molecular dynamics simulations performed 
with water [45]. However the experimental slip lengths 
show a systematic bias toward higher values than theo- 
retical ones. The discrepancy is actually not understood 
and might be the effect of gas adsorbed at the wall, pos- 
sibly condensed in nanobubbles, an active reseach area 
recently [M1IH5]. 

Finally, the case of organic molecules still raises de- 
bates. In the case of perfect wetting, various works on 
mica surfaces [T31 [TB] , as well as on glass surfaces [TJ] [7D] , 
showed that alkanes and OMCTS have a no-slip velocity 
either at the wall or at a molecular size inside the liq- 
uid phase (negative slip lengh) . However Pit et al. , using 
a ffuorescence recovery method with molecular probes 
in a TIRF configuration, found very large slip lengths 
(~ lOOnm) of hexadecane and squalane on atomically 
smooth saphir surfaces although these liquids wet per- 
fectly the saphir [57] [55] . Further work is under progress 
to understand if these large slip lengths in the case of 
favorable solid-liquid interactions refiect a particular ef- 
fect of molecular structure and incommensurability at 
the saphir/liquid interface. More generally, although the 
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problem of the boundary condition has generated an im- 
posant instrumental progress in the quantitative mea- 
surement and control of flows at a nanoscale, a full quan- 
titative confrontation with the theory, e.g. Eq.(lll, is 
still lacking. 



IV. NANOMETRIC DIFFUSE LAYERS AND 
TRANSPORT THEREIN 



The transport of electrolytes within the EDL and the 
associated phenomena have been extensively studied dur- 
ing the last century due in particular to its central role in 
colloid science. It has been accordingly discussed exhaus- 
tively in textbooks and reviews [lH [SH] and in the present 
section we merely point to the main concepts which will 
be useful for our discussion. 

In spite of this long history, new directions and per- 
spectives have been opened very recently in this old do- 
main. Novel insight into the nanometric Debye layer have 
indeed been reached: new tools now make it possible to 
investigate the detailed structure and dynamics inside 
the Debye layer. This concerns both the experimental 
side, see previous section, and numerical side, with the 
development of ever more powerful Molecular Dynamics 
tools. 

From these more 'molecular' views of the EDL, in con- 
trast to the more standard 'continuum' picture of the 
EDL, new concepts and phenomena have emerged. We 
shall thus merely focus here on the recent progress and 
new leads opened by nanofluidics in this domain. 



A. Prom Debye layer to ion specific effects 

An EDL builds up at a charged interface: it is the 
region of finite width where the surface charge is bal- 
anced by a diffuse cloud of counterions [55]. The Debye 
layer is usually described on the basis of the mean-field 
Poisson-Boltzmann (PB) theory, see [39, for a review. 
Ions are described as point charges and only their valency 
enters the description. Correlations between charges are 
accordingly neglected. The thermodynamic equilibrium 
balancing the ions entropy to their electrostatic interac- 
tion with the surface (attraction for the counter-ions and 
repulsion for the co-ions) leads to the introduction of the 
Debye length, as introduced above in Eq. (|5| and which 
we recall here: 



A, 



-1/2 



(14) 



Getting more into details, the PB description results 
from the combination of (i) the Poisson equation relat- 
ing the electrostatic potential V to the charge density 
Pe = e(p^ ^ P~) according to 



with e the dielectric constant, here identified to its bulk 
value, and A the spatial Laplacian; and (ii) the ther- 
modynamic equilibrium leading to a spatially constant 
electro-chemical potential for the ions: 



M± = Ai(p±) ± z±eV 



(16) 



Ay 



e 



(15) 



with z± the ion valency, pit the ion densities and under 
the above assumptions, /^(p) = /io + ^sTlog p the perfect 
gas expression for the chemical potential. 

Though its crude underlying assumptions, the mean- 
field PB theory captures most of the physics associated 
with the EDL and is the usual basis to interprete ion 
transport. One reason for this success is that ion cor- 
relations within the EDL, which are discarded in PB, 
can generally be neglected. This can be quantified by 
introducing a coupling parameter. Tec, which compares 
the typical inter-ions electrostatic interaction to thermal 
energy [Ml El]. Following Refs [90] and [91], the latter 
is estimated as T^c « / {Aiitl) /ksT = is/i, with i 
the mean distance between ions at the surface. This 
scale is related to the charge density E, as i w E"^/^, 
so that Fee ~ \/'^f-B El- is expected to break 

down when the coupling parameter exceeds unity. For 
typical surfaces this threshold is not reached: for glass, 
the surface charge is at most in the 10~^C/m^ range 
(^ G.IO^'^ e/nm^), and the coupling parameter F^c re- 
mains smaller than unity. 

Other assumptions inherent to PB could be also ques- 
tioned. In the PB description, the dielectric constant is 
assumed to be spatially uniform and identical to its bulk 
value. This assumption is expected to break down very 
close to surfaces, in particular hydrophobic, where water 
ordering induce a local electric dipole pointing outward 
to the interface [HI]. This casts some doubts on the local 
relation between medium polarization and local electric 
field. However MD simulations suggest that using the 
bulk dielectric constant in Poisson equation is not a crit- 
ical assumption and can be safely used to describe the 
EDL [92 . 

PB description however misses an important class of 
effects associated with the specific nature of ions, and 
which affect the fine structure of the EDL. Ion specificity 
is intimately connected to Hoffmeister effects, namely 
that the interactions between charged and neutral objects 
in aqueous media do depend crucially on the type of ion 
and not only on electrolyte concentration [33] [M] [55] . 
Such effects occurs not only at the air-water interface, 
where it does affect surface tension depending on ion 
type, but also at hydrophobic surfaces ^V, "M*. This 
effect is evidenced in Fig. [5] where results of molecular 
dynamics simulations of ions dispersed in water (SPC-E 
model) from Ref. [96] demonstrate that iodide ions do 
strongly adsorb at air-water and hydrophobic interfaces, 
while chloride do not. Note that a EDL builds up even 
in the case of a neutral interface, we shall come back on 
this result and its implications. In general heavier halide 
ions adsorb at the hydrophobic surface. This effect is 
mostly absent at a hydrophilic interface. Tools of molec- 
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FIG. 5: Simulated density profiles of negative ions (solid 
lines), positive ions (dashed lines), and water (dotted lines) 
for roughly 1-M solutions of (a) NaT and (b) NaCl between 
neutral hydrophobic surfaces, (c) Nal at a vapor-liquid inter- 
face, and (d) Nal between neutral hydrophilic surfaces. From 
Ref. [Ml- 



ular simulation and new spectroscopic techniques have 
fully renewed the interest in ion specific phenomena |95j . 
The origin of the effect involves many aspects, as ion 
polarizability, image-charge interaction, ion-size (steric) 
effects, dispersion forces and the mechanisms are still the 
object of intense research ... and also of some fierce de- 
bate, such as for the prefered adsorption of hydronium 
versus hydroxyl at interfaces |97j. which fixes the charge 
of the liquid- vapor interface. 

Such ion-specificity effects, which go beyond the tradi- 
tional PB framework, have profound effects on the ion- 
transport process within the Debye Layer, which we ex- 
plore in the next part. 



B. Interfacial transport: Electro-osmosis , 
streaming currents, and beyond 

The EDL is at the origin of various electrokinetic phe- 
nomena. Here we shall discuss more particularly electro- 
osmosis (and -phoresis for colloids) and its symetric phe- 
nomena, streaming currents. Both phenomena take their 
origin in the ion dynamics within the EDL. Due to 
the nanometric width of the EDL, these phenomena are 
^nanofluidic' by construction and the various lengths in- 
troduced in Sec will appear. 

Let us first discuss electro-osmosis (EO). EO is the 
phenomenon by which liquid flow is induced by an electric 
field. In its simplest form, the fluid velocity Vf beyond 
the EDL is linearly connected to the applied electric field 
Ef. according to the Smoluchowki formula |42j : 



(17) 



with e the dielectric constant, 77 the viscosity and C the 
so-called zeta (electrostatic) potential. The common in- 
terpretation of C is that it is the electrostatic potential 
at the "shear plane", i.e. the position close to the wall 



where the hydrodynamic velocity vanishes. As we now 
discuss, this interpretation should be taken with caution. 

Let us briefly recall the derivation of Eq{T7] Fluid dy- 
namics is described by the stationary Stokes equation: 



PeEe = 



(18) 



Here we use a frame where the electric field i?e is along x 
and parallel to the planar surface, and the velocity field 
Vx only depends on the direction z perpendicular to it. 

A key remark, already mentionned above, is that the 
charge density, = e(p^ — p_), is non vanishing only 
within the EDL, so that the driving force for fluid mo- 
tion, Fe — PeEe is limited to that nanometric region, and 
vanishes otherwise. In order to obtain the velocity pro- 
flle, the Stokes equation should be integrated twice, and 
thus requires two boundary conditions. Far from the sur- 
face, the velocity profile is plug-like and 9zv|z=oo — 0. At 
the wall surface, the hydrodynamic boundary conditions 
for the hydrodynamic velocity profile should be specified. 
We have discussed above that the latter introduces a slip 
length 6, according to Eq. [t] as bdv\z=o = v{z = 0). Us- 
ing Stokes and Poisson equations, one gets immediately 
that the full velocity profile is related to the electrostatic 
potential V as 



v{z) = --E X [-V{z) + C] 



(19) 



where C has the meaning of a zeta potential |153j and 
takes the expression 



C^VoX{l + b- Keff) 



(20) 



where Vq is the electrostatic potential at the wall and 
Kcff is the surface screening parameter, defined as k^s — 
-V'{0)/Vo [SHIM]- For weak potentials this reduces to 
the inverse Debye length: Kefr ~ A^^, but its full expres- 
sion takes into account possible non-linear effects. The 
fluid velocity far from the surface (i.e. outside the EDL) 
has the Smoluchowski expression above, Eq. (17), with 
the expression of the zeta potential, (, given in Eq. 20 
This prediction was first discussed in Ref. [99], and red- 
erived independently in more recent works fM) [55] . We 



emphasize that the above result, Eqs. (19 20 1, does not 



make any assumption on the specific model for the EDL, 
except for a constant dielectric constant e and the valid- 
ity of Stokes equation, which was discussed above to be 
valid down the nanometer. 

The physical meaning of the above expressions is quite 
clear. The velocity in the fluid results from a balance 
between the driving electric force and the viscous friction 
force on the surface. Per unit surface, this balance writes 



rj 



eE X E 



(21) 



with S the surface charge density (in units of rn^'^) at 
the wall. The expression for the viscous friction expresses 
the fact that flow gradients occur on a length b + Ad . 
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FIG. 6: Sketch of the influence of shppage on the electro- 
osmotic transport. Slippage reduces the viscous friction in the 
electric Debye layer, as the hydrodynamic velocity gradient 
occurs on a length b + Xd, instead of A without slippage. 
Flow is accordingly enhanced by a factor 1 + b/Xo- 



Relating the surface charge to the surface potential Vq 
(eE = —edzVlo), one recovers the above results. 

As a side remark, one may note that the slip enhance- 
ment does not require any assumption on the electric and 
fluid behavior. Indeed at a fully general level, the (ex- 
act) momentum balance at the surface shows that the 
velocity at the wall is exactly given by Uwaii = -^p" E, 
with I?waii the dielectric displacement computed at the 
surface, and A the surface friction coefficient introduced 
in Eq. ^ . The latter is related to slip length according to 
the definition A = b/rj. This expression for the wall veloc- 
ity under an applied electric field neither invokes Navier- 
Stokes equation, nor any assumption on the dielectric 
behavior. This wall velocity acts as a lower bound to the 
asymptotic EO velocity, i.e. the one measured at infinite 
distance from the surface. This demonstrates that the 
enhancement of the EO velocity due to low liquid-solid 
friction, i.e. large slippage, has a fully general validity. 

Several conclusions emerge from the above equation. 
First if the wall is characterized by a no-slip BC, with 
b — 0, then the zeta potential Eq- ( pOj ) identifies with 
the electrostatic potential Vq computed at the surface (up 
to a possible molecular shift due to the precise position of 
the hydrodynamic BC which defines the so-called shear 
plane [98]): this is the common assumption, as usually 
expressed in textbooks [12]. Now, if finite slippage oc- 
curs at the surface, then the zeta potential is much larger 
than the surface electrostatic potential, Vq: it is amplified 
by a factor 1 + b/Xjj (assuming for the discussion that 
Kcff = A^^). This factor may be very large, since on bare 
hydrophobic surfaces b may reach a few ten nanometers 
(say b ^ 20 ~ 30nm), while the Debye lentgh typically 
ranges between 30 nm down to 0.3 nm (for 10~'*M to 
IM). 

These predictions have been confirmed by molecular 
Dynamics simulations of electro-osmosis (and streaming 
current, see below) [92, 96, 98}. Various systems were 
considered, from a model electrolyte involving ions in a 
Lennard- Jones solvent, to a more sophisticated SPC-E 
model of water, with ions of various nature. Altogether 
simulations confirm the key role of slippage on electroki- 
netics and results do fully confirm the above description. 



On the experimental side, this strong influence of sur- 
face dynamics on the electrokinetics at surfaces has not 
been appreciated and explored in the electrokinetic lit- 
terature up to now, and very few experimental investi- 
gations have been performed on this question [ 1001 llOl] . 
One key difficulty is that the above result for the zeta po- 
tential, Eq. (20 1, involves a strong entanglement between 
electrostatics, through surface potential, and fluid dy- 
namics, through hydrodynamic slippage. In order to dis- 
entangle the two effects, two independent measurements 
of these quantities should in principle be performed. To 
our knowledge the first work on the subject was per- 
formed by Churaev et al. in Ref. |100j . and results 
indeed suggested a slip effect, at the expense however 
of a rather uncontrolled assumption on the surface po- 
tential. More recently this problem was tackled in Ref. 
|101j using the nanoPIV tool discussed above. This set- 
up allows for two independent measurements of C and 
Vq, leading to an unambiguous confirmation of the above 
predictions of slippage effect on the zeta potential. In- 



deed both the velocity profile, which is fitted to Eq. 19 
and the nanocoUoid concentration profile, from which the 
surface potential is deduced on the basis of the electro- 
static repulsion, are measured independently. Velocity 
profiles from Ref. |101j are displayed in Fig. [t] While the 
two surfaces under consideration (a hydrophilic glass and 
a silanized hydrophobic surface) had basically the same 
surface potential, a factor of two is found on the zeta po- 
tential under the conditions of the experiment. Results 
are consistent with a slip length of w 40nm. The fac- 
tor of two for C/^o occurs here due to a Debye length 
of ~ 50nm in the experimental condition of Ref. |101j 
(a large Debye length is required to investigate fiow in- 
side the Debye layer at the present spatial resolution, 
^ 20nm, of the nanoPIV technique). As a side remark, 
note the large slip velocity at the wall in the slippy case 
in FigjT] this is expected since according to EqjT9] the 
velocity at the surface is directly proportional to the slip 
length b, as v{0) ~ ^fj^ ^ ^ (assuming, to simplify, 

that Kcfi ~ A^^"'^). 

This result opens new perspective to make use of 
surface physico-chemistry in order to optimize electric- 
induced transport and control flow by surface properties. 
It was also argued that such slip-induced optimization 
would strongly enhance the efficiency of energy conver- 
sion devices based on electrokinetic effects |1021 1103j . 

Besides slip effect, ion-specificity has also an impor- 
tant influence on the electro-osmotic transport. As we 
reported above, the structure of the EDL is affected by 
the nature of the ions under consideration, especially at 
hydrophobic walls (and air- water interfaces). This can 
be best seen by rewriting the expression for the zeta po- 
tential, in Eq|20] in a slightly different form: 



c 



1 



{z' + b)pe{z')dz' 



(22) 



with Pe the charge density. This expression e.g. comes 
from the double integration of the Stokes law. As seen in 
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FIG. 7: NanoPIV measurements for the electroosmotic flows 
close to hydrophilic pyrex • and hydrophobic OTS-coated sur- 
face ■ . The deduced zeta potentials are C, — — 66± 8mV (hy- 
drophilic glass) and ^ = — 123±15mV (hydrophobic, silanized 
surface), while the electrostatic potentials were independently 
measured to be comparable: Vo — — 69mV for the hydrophilic 
(pyrex) surface, and Vo ~ — 65mV for the hydrophobic (OTS- 
coated) surface. The salt concentration is 1.5 10~^mol.L~'^ 
and the driving electric field is -E = 500V.m^^. Dashed lines 
are fits to the theoretical predictions using Eq. (19 1, with 
Xn = 51±10 nm (measured independently) and & = 0± lOnm 
(bottom) and 6 = 38 ± 6nm (top) for the slip lengths. Figure 
from Ref. [TOT] . 



Fig. [5] ion specific effects have a strong influence on the 
charge distribution, and do indeed strongly modify the 
value of This effect was studied in various MD simula- 
tions of electrokinetics [92l |96l I104j and several counter- 
intuitive effects were indeed observed, such as flow rever- 
sal (as compared to the expected surface charge) |104j . or 
even the existence of a non-vanishing zeta potential for 
a neutral surface [92j [96] ! These effects can be under- 
stood from Eq. [22] while for a vanishing surface charge, 
the system is electro-neutral and Pe{z') dz' — 0, this 
however does not imply that its first spatial moment, 
which enters the expression of C according to Eq. [22] be 
zero. 

More quantitatively, effects of ion-specificity on the 
zeta potential were rationalized in Ref. [ST on a basis of 
a minimal model, which was subsequently validated on 
MD simulations. This description is based on the idea 
that 'big' ions have a larger solvation energy and are 
therefore attracted to hydrophobic interfaces in order to 
minimize this cost |105j . An estimate of the solvation 
free-energy of ions close to interface allows to propose 
analytical predictions for ion-specific effects on the zeta 
potential. 

We finally quote that all the above discussion on the 
electro-osmotic transport also applies to the streaming 
current phenomenon, by which an electric current is in- 
duced by the application of a pressure gradient. Again, 
since the EDL is not electro-neutral, liquid flow will in- 
duce an electric current inside the EDL [42 . According 



to the Onsager symetry principle, these two phenomena 
- electro-osmosis and streaming current - are intimately 
related and the cross coefficients should be identical. The 
expression for the streaming current under an applied 
pressure gradient VP takes the form: 



Ie = ~—A{-VP) 



(23) 



where C takes the same expression as in Eq.(20l and A 
is the cross area of the channel. 

We end up this section by mentioning that the above 
effects generalize to any interfacial transport phenom- 
ena. Indeed the above electrically induced phenomena 
belong to a more general class of surface induced trans- 
port, which also involves phenomena known as diffusio- 
osmosis or thermo-osmosis and their associated phoretic 
phenomena for colloidal transport [2TJ |89]. Thermo- 
osmosis points to fluid motion induced by thermal gra- 
dients, while diffusio-omosis corresponds to fluid flow in- 
duced by gradients in a solute concentration. Though 
relatively old, thermophoretic transport has been the ob- 
ject of recent investigations in particular for colloid ma- 
nipulation |106l 11071 1108] . Its origin is still the object 
of an intense debate with recent progress in disentan- 
gling the various contributions |106[ 11071 1109j . On the 
other hand, diffusiophoretic and osmotic transport have 
been less explored [SHI IllOj but their potential in the con- 
text of microfluidic applications has been demonstrated 
recently |11H 1112] . This opens novel possibility of driv- 
ing and pumping fluids |112j - as well as manipulating 
colloids - with solute contrasts, which have been barely 
explored up to now and would deserve further investiga- 
tions in the context of nanofluidics. A common point to 
all these phenomena is their "nanofluidic root" : as for 
the electro-osmotic transport discussed above, the driv- 
ing force for fluid motion is localized within a diffuse layer 
of nanometric size close to the surfaces. Therefore, this 
opens the possibility of strongly amplifying their effect 
on the basis of slippage effects at the solid interface. As 
demonstrated theoretically in Ref. |113j . the amplifica- 
tion is - as for electro-osmosis above - amplified by a 
factor 1 + b/X, where A is the width of the diffuse layer. 
Similar effects are predicted for thermophoretic transport 

Finally, one may raise the question of interfacial 
transport on super-hydrophobic surfaces. Such surfaces, 
achieved using nano- or micro- engineering of the sur- 
faces, were shown to considerably enhance the slippage 
effect and exhibit very large slip length in the micron 
(or even larger) range j46t 1471 UlSj . Accordingly, a naive 
application of the previous ideas would suggest massive 
amplification by a factor up to 10* (!). However, the 
composite structure of the superhydrophobic interface, 
involving both solid-gas and liquid-gas interfaces, makes 
these transport mechanisms far more complex in this case 
than on smooth interface. It has been shown that in 
the regime of thin Debye layer, no amplification is ob- 
tained for electro-osmosis on superhydrophobic surfaces 
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|112l 1145] . unless a non- vanishing charge exists at the 
liquid-gas interface. In contrast, a massive amplification 
is predicted for difi^usio-osmosis |112j . This strong pre- 
diction has not received an experimental confirmation up 
to now. 

Furthermore, a massive amplification, by a factor up 
to 10^ (!), may be achieved when using nano-engineered 
surfaces associated with super-hydrophobic properties, 
which considerably enhance the slippage effect (with b 
raising to the micrometers range) [i51H71[ll5] . 



C. Surface versus bulk: conductance effects 

We now discuss the surface conductance effect in 
nanochannels. 

The conductance, K, characterizes by definition the 
electric current versus electric potential drop relation- 
ship. As pointed out previously, conductance probes the 
number of free charge carriers, here ions. In the bulk 
K is therefore expected to be proportional to the salt 
concentration ps (ionic strength). Now in the presence 
of surfaces, the charges brought by the surface lead to a 
supplementary contribution to the conductance. 

Surface conductance is actually a rather classical phe- 
nomenon in colloid science [12 , but it was demonstrated 
quite recently in the context of nanochannel transport by 
Stein et al. [IT] and then by other groups |H] . 

As we introduced in Sec. ]IIB 2| c, the ratio of bulk 
to surface charge carriers is characterized by the Dukhin 
length £du- This length accordingly describes the com- 
petition between bulk and surface conductance. Depend- 
ing on salt concentration, this length can be much larger 
than the Debye length. Surface conductance is there- 
fore expected to dominate over the bulk contribution for 
nanochannels smaller than £du- The magnitude of the 
surface conductance effect is amplified in small channels, 
but a key point is that it does not require Debye layer 
overlap. 

Experimentally, the surface contribution to the con- 
ductance shows up as a saturation of the conductance in 
the limit of small salt concentration, while its expected 
bulk counterpart is expected to vanish in this limit. The 
saturation originates in the charge carriers brought by 
the surface charge, the number of which is independent 
of the salt concentration. This is illustrated in Fig. [8] 
from [41 , where the conductance is measured in chan- 
nels of various width, from 1015nm down to 70nm: 

Let us discuss more specifically this point. In a slit 
geometry with width h, the general expression for the 
current (for a unit depth of the slit) takes the form 

rh/2 

Ie = e dz [p+(z)u+(z) - /3_(z)u_(z)] (24) 



where the averaged velocites of the ions is 
u±{z) = v{z) ± efi±Ee, 
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FIG. 8: Channel height dependence of ionic conductance be- 
havior. The conductance of fluidic channels is plotted against 
salt concentration for various nanochannel width h {h = 1015, 
590, 380, 180, and 70 nm). The inset displays the fit values 
for the surface charge at the nanochannel walls as a function 
of h, from Ref. E]. 



with fj,± is the ion mobility, defined here as the inverse 
of the ion friction coefficient (we assume furthermore 
= m)- The velocity v{z) is the fluid velocity induced 
by electro-osmosis under the electric field, as given in 
Eq.(|l9|). The second term is the contribution to the cur- 
rent due to ion electrophoresis. 

Combining the expression of the velocity v{z), as 
( 24p5 ), one may remark that the 



in Eq.(19l, to Eqs 



electro-osmotic contribution to the current takes the form 
of the integral of the charge density times the electro- 
static potential: this has therefore the form of an electro- 
static energy, £. This is confirmed by the full calcula- 
tion, which leads to a conductance expression in the form 

mm-- 



(26) 



with A the channel cross section (A — w x h with ?« the 
depth of the channel) and the surface correction H 



H = cosh 



eVc 



£ 



ipsksTh 



1 



(27) 



where Vc — V{0) is the electrostatic potential at the cen- 
ter of the channel, and £ — ^ /siit '^'^ (iz")^ electro- 
static energy (and not the free energy). Note that writ- 
ing p = Snrjcli with di the ion diameter, the last term in 
Eq.(27) is typically of order di/£B- It originates in the 
surface-induced - electro-osmotic - contribution to the 
conductance. 

In the low salt concentration regime, p^ — )■ 0, the 
length scales are organized in the order: 



^GC < Xd < (-Du, 



(28) 



(25) 



since XdI^gc = (-Du/^d S> 1 for pg 0. Saturation will 
occur as the channel width h is smaller than the Dukhin 
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length, h < iou, independently of the order of the h and 
Xd- Debye layer overlap is accordingly not the source of 
the saturation of the conductance. In this regime, the 
second term in Eq. (27 1, associated with surface contri- 
butions is dominant. 

In this limit, the calculation of the electrostatic energy 
shows that £ ~ The fact that £ scales linearly in 
and not as |Sp - as would be first guessed - is due to the 
non-linear contributions to £ in the Poisson-Boltzmann 
description. These become dominant in the ps limit 
where the Debye length becomes larger than the Gouy- 
Chapman length. Altogether, this shows that for ps 0, 
the conductance saturates at a value 



pw X 2|S| 1 



1 



27r^s p rj 



(29) 



with w the depth of the channel. This value is inde- 
pendent of both h and Ps- This saturation can be also 
partly understood in the context of Donnan equilibrium 
discussed below. In the limit ps 0, only the counter- 
ion contributes to the ion concentration in the channel, 
-I- p_ K, 2S//i, see Eq.(34l below, so that the con- 
ductance reduces to Kgat ~ e^^w2|S|. This analysis 
however misses the electro-osmotic contribution to the 



conductance (second term in the brackets in Eq. (29 1). 

A few remarks are in order. First, the dependence of 
the conductance on the surface properties, here the sur- 
face charge S, opens new strategies to tune the nanochan- 
nel transport properties, via an external control. This 
has been used by Karnik et al. to develop a nanofluidic 
transistor [5], allowing to control the conductance of the 
nanochannel thanks to a gate voltage. 

Another interesting remark is that along the same line 
as in Sec. |IV B| the electro-osmotic contribution to the 
conductance could be also amplified by slippage effects. 
This can be readily demonstrated by a direct integra- 
tion of Eq. (p4| . Accordingly slippage effects add a new 



contribution to the conductance: 



Ks\m = 2e w 



X b 



V 



(30) 



Assuming Stokes law for ion inverse mobility {p — 
(STrrydion)^^ with dion the ion size), then /Cgiip/^sat ~ 
b X dion |S| « b/^GC- since the Gouy-Chapman length 
£gc is typically nanometric (or less), this enhancement 
is therefore very large, even for moderate slip length 
(6 ~ aOnm) ! 

The channel resistance is accordingly decreased by the 
same factor in the low salt regime. This amplification of 
conductance by hydrodynamic slippage opens very inter- 
esting perspectives in the context of electrokinetic energy 
conversion, as was pointed out recently by Pennathur et 
al. |103j and Ren and Stein [102j. Moderate slippage, 
with slip length of a few ten nanometers, is predicted to 
increase the efficiency of the energy conversion up to 40% 
(and of course even more with larger slip lengths). This 
attracting result would deserve a thorough experimental 
confirmation. 



D. Debye layer overlap and nanofluidic transport 

The phenomenon of Debye layer overlap has already 
been widely explored in the nanofluidic litterature. As 
we pointed above, various reasons underlie this specific 
interest: (i) well controlled nanometric pores with size 
in the range of 20-100 nm can be produced using micro- 
lithography techniques |40j : this does indeed correspond 
precisely to the range of typical Debye lengths for usual 
salt concentrations (remember that Xd = 30nm for a salt 
concentration of lO^^'M); (ii) the overlap of Debye length 
does indeed have a strong influence on ion transport, so 
that novel transport effects emerge at this scale, with 
applications for chemical analysis. 

The question of Debye layer overlap, and related phe- 
nomena, was discussed quite extensively in a recent re- 
view by Schoch et al. [40]. Here we thus focus on the 
main guiding lines underlying this problem and point to 
the key physical phenomena at play. In particular we 
shall discuss an illuminating analogy to transport in semi- 
conductors. 



1. Donnan equilibrium 

An important notion underlying the Debye overlap is 
the so-called "Donnan equilibrium" , which is a well know 
concept in the colloid litterature. Due to the supple- 
mentary charges brought by the nanochannel's surface, a 
potential drop builds up between the nanochannel's in- 
terior and the external reservoir, in order to maintain 
a spatially uniform chemical potential of the ions. The 
latter, as introduced above, takes the form (assuming a 
dilute ion system): 

Kp±) = Po + kBTlogip±)±eV = po + ksTlogips) (31) 

with Ps the (uniform) salt concentration outside the 
nanochannel, i.e. the "reservoir". These equations are 
completed with the overall electroneutrality over the 
channel, 



h/2 



dz (p+(z)-p_(z)) = 2E 



(32) 



-h/2 



with z along the direction perpendicular to the nanochan- 
nel (here assumed to be a slit), see Fig. |9] 

We now focus on the thin pore limit, assuming Debye 
layer overlap h < 2Xd. In this situation, ion densities 
and electrostatic potential are approximately spatially 
uniform over the pore thickness h. Ion densities thus 
obey: 



p+x p_= p^ 
2S 

P+- P- ^ -r 



(33) 
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FIG. 9: (a): Sketch of the electrostatic potential in a charged 
slit with thickness h: (left) h > 2\ij; and (right), h < 
2Ad, corresponding to Debye layer overlap, (b) A charged 
nanochannel, with thickness h and length L, connects two 
reservoirs Ri and R2- In the situation of Debye layer overlap, 
a Donnan potential Vd builds up along the nanochannel. 



Altogether, this leads to 



P± 




(34) 



where the mean electrostatic potential Vd is denoted as 
the Donnan potential. 

Note that the amplitude of the Donnan potential Vd 
is quantified by the ratio iDu/h where Idu = "^/Ps is the 
Dukhin length introduced above. To fix ideas Idu ^ Inm 
for ps — IM and a typical surface charge. 

It is finally of interest to mention the strong analogy 
between the Donnan equilibrium and the equilibrium of 
charge carriers in doped semi-conductors (SC). In the 
latter the electrons density /?„ and the holes density Pp 
obey Eq. (33 1 with the electrolyte density ps replaced 



by the carriers density in the intrinsic (non-doped) SC. 
The surface charge in the nano-channel actually acts as 
impurities delivering additionnal carriers, with S > 
corresponding to donnors (N-doped SC) and S < ac- 
ceptors (P-doped SC) . The Donnan potential is then the 
analogous of the shift in the Fermi energy due to impu- 
rities, which determines the voltage difference at equilib- 
rium between regions of different dopage level. The table 
hereafter, table I, summarizes the corresponding quanti- 
ties in the nano-channel and SC analogy. 

The Donnan potential and Debye layer overlap have 
many implications on nanofluidic transport. 



Nanofluidics 


Semi-conductors 


P— 

negative ions concentration 


Pn 

electron density in the conduction band 


P+ 

positive ions concentration 


Pp 

holes density in the conduction band 


Ps 

electrolyte concentration 


carriers density in the intrinsic SC 


surface charge 211 /h 


Impurities (dopping) concentration 


Donnan potential 


Shift of the Fermi level 



TABLE I: Equivalence table for corresponding quantities in 
doped semiconductors and nano-channels. 



2. Ion transport and PNP equations 

A classical framework to discuss ion transport in nar- 
row pores is the so-called Poisson-Nernst-Planck (PNP) 
equations. This simplified description of ion transport 
in strongly confined channels is based on the coupled 
diffusion-electro-convection of the ions |117j . 

As in Figj9] the channel thickness is assumed to be 
small compared to the Debye length, h <C A^, so that 
ion concentrations are assumed to be uniform across its 
thickness. The ion fluxes J± (per unit surface) then takes 
the form 

J± = -D±dxp±{x) ± p±e p±{x){-d^V){x) (35) 

where D± and p,± are related by the Einstein relation- 
ship, D± = fi±kBT; X is along the direction of the chan- 
nel. Note that for simplicity, a single salt and mono- 
valent ions are considered here. A local version of the 
electroneutrality, Eq. (33 1, may be made along the pore: 



2E(x) 
h 



(36) 



Finally at the pore entrance and exit, a Donnan elec- 
tric potential drop builds up, along the description given 
above, see Eq. (34 1. 

To illustrate further this approach, we consider the 
transport through a single nanochannel, as sketched in 
Fig|9]-b. At equilibrium the two reservoirs, i?i and R2 
have the same electric potential and salt concentration, 
associated with a uniform Donnan potential along the 
channel. Now if an electrostatic potential drop AT^ or a 
salt concentration difference is imposed between the 
reservoirs, ion fluxes will build-up so as to relax towards 
equilibrium. 

In the stationary state, the ion flux are spatially uni- 
form, so that one may solve Eqs. ( 35]|36 1 for the ion 
densities and electrostatic potential. This allows to com- 
pute the current / — Ae{Jj^ — J_) and total ion flux 
<i>t — A{J+ + J-), with A the cross area. In general, the 
relationship /(AF, Aps), <I>t(Ay, Ap^) are non-linear. 

However, in the limit of small AV^ and Aps, a linear 
relationship can be written in general between the ionic 
fluxes and corresponding thermodynamic 'forces' |119l 



(37) 



" I ' 


_ A 


■ K 


P'K 




-AT/ 




^ L 




Moff 




-fcBTA[logp, 



18 



with A the cross area of the channel and A[logps] = 
Aps/ps- Note that due to Onsager (time reversal) 
symetry, the non diagonal coefficients of the matrix are 
equal [TTmiT^ . 

The coefficients of the above matrix can be calculated 
within the PNP framework. PNP equations are solved 
analytically for small potential and concentration drops, 
leading to the following expression for the various coeffi- 
cients in the matrix, Eq.(37l: 



^off = 2p.x 



p.K ^ ep.x 




(38) 



We assumed here that ions have the same mobilites fi± = 
II. The cross effects, associated with the mobility p,K, 
originate in the dependence of the Donnan potential on 
the salt concentration. 



Let us conclude this part with a few remarks: 

• First one should realize that a number of assump- 
tions are implicitly made in writing PNP equations. In 
particular the ions are treated as a perfect gas and cor- 
relations between ions along the channel are neglected. 
This may become problematic in strongly confined 
situations where single file transport (of the solvent) and 
strong unidirectional electrostatic correlations should 
build up. As an example, proton transport in single file 
water has been shown to involve a highly-cooperative 
mechanism |118j . However these limitations are re- 
stricted to single file transport in molecular channels, 
and should not be a limitation for pore size larger 
than a nanometer. Furthermore the PNP model is 
interesting to explore as guiding line, in the sense that it 
provides a rather correct physical idea of the (complex) 
electro-diffusion couplings. 

• Furthermore, we may pursue the analogy discussed 
above between nanochannels and doped SC (see Ta- 
ble I), and extend the discussion to transport phe- 
nomena. Indeed the Poisson-Nernst-Planck equa- 
tions, Eqs.( [35| ), are formally identical to the phe- 
nomenological transport equations for electrons 
and holes in SC. Thus, as long as electrical and concen- 
tration fluxes only are allowed, nanofluidic devices can in 
principle reproduce standard SC-based components, such 
as diodes and transistors. The nanofluidic diode for in- 
stance is based on the properties of a PN junction, i.e. 
the junction between two regions with different doping: 
according to the equivalence table, Tabel I, its nanofluidic 
equivalent corresponds to two nanochannels with differ- 
ent surface charge S/ft. We will discuss below (paragraph 



5) some recent findings confirming the pertinence of this 
analogy. 

• However one should keep in mind that the trans- 
port analogy between fluids and electrons breaks down 
in the presence of an hydrodynamic flow. Actually, PNP 
equations introduced above do not take into acccount 
convective contributions. In a fluidic system, a flow is 
indeed expected to occur in nanochannels as soon as a 
voltage drop is applied to its end, due to the body force 
acting on the mobile charges, but also under a pressure 
gradient (hydrostatic or osmotic). This is for example 
the source of electro-osmotic or streaming current con- 
tributions discussed above. In contrast, SCs can not flow 

- quite obviously - under an applied electric field and may 
only deform elastically. This is therefore a limitation to 
the above analogy. 

Beyond this analogy to SC transport, the above PNP 
equations do not involve any hydrodynamic contribu- 
tions. This concerns for example the electro-osmotic con- 
tributions (see e.g. the very last - viscosity dependent 

- term in the expression of the conductance, Eq.(27)), 
or the slip-induced enhancement of these phenomena. 
These may well become dominant in the strongly con- 
fined regime for sufficiently large slippage. In general 
a convection contribution to the ion flux should there- 
fore be added to the above ion flux, in the form J^^v ~ 
VfXp±, where Vf is the water velocity. The water velocity 
is in turn coupled to (electric) body forces via the mo- 
mentum conservation equations, in practice the Navier- 
Stokes equation in its domain of validity. Altogether, the 
hydrodynamically induced cross phenomena, i.e. electro- 
osmosis, streaming currents, diffusio-osmosis, ... (as well 
as bare osmosis) , can be summarized at the linear level in 
terms of a symetric transport matrix, in a form similar to 
Eq. ( 37 1 . This matrix relates linearly the fluxes to ther- 
modynamic forces at work |120l 11211 I122j . These phe- 



nomena have been discussed extensively in the context 
of membrane transport but generalize here to transport 
in nanofluidic channels. 

More generally the coupling of the ion transport to the 
fluid transport thus opens the route to new applications 
in nanofluidic devices, with a richer phenomenology than 
in SC electronics. 

We now illustrate the above concepts on a few practical 
situations. 



3. Perms electivity 

As a first example, we consider the permeability of 
nanopore to ions. It was first shown by Plecis et al. 
[2] that nanochannels exhibit a selective permeability 
for ion diffusive transport, see Fig. 10 'a): ions of the 
same charge as the nanochannel surface (co-ions) exhibit 
a lower permeability, while ions of the opposite charge 
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(counter-ions) have a higher permeability through the 
nanochannel. 

This charge specific transport is a direct consequence 
of a non-vanishing Donnan potential in the nanochan- 
nel (and Debyc layer overlap). As pointed out above, 
Fig. [9] counter ions exhibit a higher concentration in 
the nanochannel as compared to that of a neutral specie, 
/9+ > Pn, while co-ions have correspondingly a lower con- 
centration, /9_ < Pn- 

The diffusive flux of counter-ion will be accordingly 
larger than its expectation for neutral species, and vice- 
versa for co-ions. This leads therefore to a charge-specific 
effective diffusion coefficient Dcs for the co- and counter- 



with p±{x) the (local) ion concentration in the nano- 
channel, averaged over the channel width h; q the ion 
charge; p{x) the local salt concentration in absence of 
electrostatic interactions (fixed by the salt concentra- 
tion in the reservoirs); and V{z) the electrostatic po- 
tential across the channel, for which Plecis et al. used a 
(linearized) Poisson-Boltzmann expression. A very good 
agreement with experimental results is found, as shown 



in Fig. 10 showing that this exclusion-enrichement pic- 
ture does capture the essential ingredients of the ion dif 
fusive transport. Note that using the Donnan descrip- 
tion above, one may furthermore approximate the elec- 
trostatic potential by its Donnan expression in Eq.(34l 
This leads to an analytical expression for /3± as 
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■ experimental neutral 
• experimental dianion 
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p(x) 




p± 
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(41) 



This 



where, as shown in Eq.(34l, p± = yp^ + (ff)^ 
expression reproduces the results in Fig. |40] 

Finally, one should note however that due to the differ- 
ence in ion permeability, a charge separation will build up 
between the two ends of the channel. This leads therefore 
to the creation of a reacting electric field along the chan- 
nel, which will compensate dynamically for this charge 
separation [2j. This points to the complex couplings as- 
sociated with charge transport in nanochannels. More 
experimental and theoretical work is certainly in order 
to get further insight in these phenomena. 



4-. Pre- concentration 



FIG. 10; Top Variation of the relative permeability of vari- 
ous probes with different charges, versus ionic strength. The 
experimental results are compared to theoretical fittings ob- 
tained with Eq. ( 40 1 (dotted lines) . Bottom Concentration 
profiles of counter- and co-ions, p±, for a given salt concen- 
tration drop A/9 between the two ends of the nanochannel. 
The dotted line represents the linear concentration profile for 
noninteracting diffusing species, whereas plain lines show the 
ion profiles in the nanochannel. Adapted from Ref. [2]. 

Following Plecis et al. [2'^, the effective diffusion coef- 
ficient is defined according to the identity 



J± 



-D 



Ap± 



-D. 



off- 



Ap 



(39) 



L ~ L 
taking into account the fact that the local concentra- 
tion in the nanochannel differs from the one imposed at 
the two ends of the reservoir. The linearized Poisson- 
Boltzmann equation is then used to calculate the ion 
concentration in the slit, p±, as: 



/3± = 
D,s/D : 



p{x) 



dz exp[TeV{z)/kBT] 



slit 



(40) 



Along the same lines, another consequence of De- 
bye layer overlap is the 'pre-concentration' phenomenon, 
which was first observed by Pu et al. [3] , and others |123j . 
Instead of a concentration drop as above, a voltage drop 
is applied along the nanochannel. It is then observed 
that ions enrich at one end and deplete at the other end. 
This phenomenon is illustrated in Fig. [TT] 

As above, a key point underlying the phenomenon is 
that ions of the opposite sign to the channel's surface, 
counter-ions, are more heavily transported than co-ions. 
This will create an enriched/depleted zone for both ions 
at the two ends of the channel. 

This phenomenon has attracted a lot of interest, due 
in particular to its potential applications in the context 
of chemical analysis, for which it would provide a very 
interesing way of enhancing the sensitivity of detection 
methods. 

At a basic level, the origin of preconcentration is a 
"cross effect" as introduced in Eq. (37 1: an electric po- 



tential drop AV leads to a ffux of ions(<i>t) , as quantified 
by the cross- mobility pK- According to the PNP result 
for pk [pk = e px2Y./h), this cross effect is thus a direct 
consequence of confinement and the existence of surface 
charges. Accordingly the magnitude of its effect, as com- 
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(a) (b) (c) (d) 



FIG. 11: Two micro-channels filled with fluorescent probes 
are connnected by a nanochannel (a). After an electric field 
is applied across the nanochannels, an ion-enrichment and 
ion-depletion occurs at each end of the nanochannel (b-e). 
Figures (a) to (e) show the evolution of the fluorescence in 
the microchannel versus time (from to 20 s). The probe 
used was (A) fluoresceine, (B) rhodamine 6G. From Ref. [3]. 



pared e.g. to the diffusive flux, depends on the Dukhin 
number, see Eq.(38). 



But the detailed mechanims underlying the preconcen- 
tration involve non-linear couplings between the various 
transport process of ions in the nanochannel, which are 
quite complex to rationalize. A systematic description 
of the phenomenon was proposed by Plecis et al. in 
Ref. f 124], revealing the existence of various preconcen- 
tration regimes. We refer to Ref. [ID] for a detailed dis- 
cussion on the mechanisms underlying this process. 



5. More complex functionalities: nanofluidic diodes 

The analogy with transport in semiconductors quoted 
above suggest that more complex nanofluid transport 
phenomena can be obtained in the regime of Debye layer 
overlap. 

We discuss here the analog of a PN junction in SC 
transport. Using the equivalence table. Table I, a 
PN junction corresponds for nanofluidic transport to a 
nanochannel exhibiting a disymetric surface charge along 
its surface. 

Such a nanofluidic diode device has been developped 
by Karnik et al. [6], following a previous work by Siwy 
et al. in a different pore geometry [7] . This is illustrated 
in Fig 12 from Ref. [5]. The surface of a nanochannel is 



coated with two different surface treatments (half with 
avidin, half with biotin), leading to a surface charge con- 
trast along the two moieties of the channel. Fig. [12]-Lefl. 
Accordingly the current versus applied electric poten- 
tial drop characteristics is found to exhibit an asymetric 
shape: as in a classical diode, the current passes only in 
one direction. Similar asymetric I-V curves are obtained 
for pores with asymetric geometries, like conical pores 
obtained by track-etch techniques [7]. 

At a more quantitative level, Karnik et al. discussed 
the effect on the basis of the PNP equations discussed 
above, under the assumption of local electroneutrality. 
In the present disymetric case, Karnik et al. solve nu- 
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FIG. 12: Left (a) A nanofluidic diode is fabricated by pattern- 
ing a nanochannel with different coatind on its two moieties 
(avidin and biotin here), (b) Epifluorescence image of the fab- 
ricated nanofluidic diode, showing fluorescently labeled avidin 
in half the channel. Scale bar 20 pm. Right current versus 
applied voltage. Symbols are experimental data, while the 
solid lines are theoretical predictions using a PNP transport 
model. From Ref. ^Gj. 



merically these PNP equations for an applied potential 
drop at the ends of the pore [B]. Solutions are shown 
in Fig. 12 demonstrating that the PNP framework is in- 
deed able to capture the transport rectification measured 
experimentally. 

This diode behavior can also be discussed in the con- 
text of the analogy with semiconductor transport that we 
put forward above. The geometry described in Fig. [12] 
is indeed equivalent to that of a classical PN junction, 
corresponding to two regions with different doping. In 
the present case of a nanofluidic diode, the doping con- 
trast between the P and N regions is associated with the 
different surface charge between the two moieties of the 
nanochannel. The diode effect in PN junctions is classi- 
cally interpreted in terms of the Space-Charge Zone at 
the interface between the two regions with different dop- 
ing. An approximate solution of the PNP equations can 
be proposed and leads to a rectified I-V characteristic 
which takes the form of the so-called Shockley equation 



S exp 



eV 
kef 



- 1 



(42) 



where Is is the so-called the saturation current Is. This 
expression is expected to describe well the blocking to 
non-blocking transition of the junction, i.e. for moderate 
eV/ksT . However this expression strongly overestimates 
the current at large eV/ksT for which a full solution of 
the PNP equations is needed. 



In the Shockley equation, Eq.(42l, the saturation cur 



rent is expressed in in terms of the donor-acceptor den- 
sities, PA and p_D, as 



Is w 



Dr, 



LnPA LpPD 



(43) 
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where rii is the carriers density, I?„, Dp the diffusion 
coefficients of electrons and holes, and L„, Lp the length 
of the P and N regions. Using the equivalence table. 
Table 1, one may convert this expression for the ionic 
transport using Z)„ = £)_, Dp — D^, Ln = Lp = L/2, 
while Pa — 2|Eyi|//i and po ~ 2\'S£i\/h. 

This Shockley expression predicts an 1-V curve in qual- 
itative agreement with the experimental result in Fig. 12 
However, quantitatively it does (very strongly) overesti- 
mate the current for large voltages V, as expected due to 
the simplifying assumptions underlying the expression in 
Eq. (|42l). 



The analogy is however interesting to capture the un- 
derlying physics behing ion rectification in ionic trans- 
port. It allows to predict a full zoology of fluidic func- 
tionalities, in line with their SC analogues. 

Finally it is interesting to quote that the disymetry 
measured here for I-V characteristics generalizes to other 
transport phenomena. For example, Siwy et al. also 
demonstrated the occurence of a disymetric diffusion of 
ions through conical pores |126l . 



Altogether these phenomena demonstrate the richness 
and complexity of behaviors obtained in the regime of 
Debye layer overlap. This definitely shows the great po- 
tential of nanofluidics, where further phenomena should 
emerge in the future taking benefit of these couplings. 



V. THERMAL FLUCTUATIONS 

In this section we raise the question of the role of ther- 
mal fluctuations in nanofluidic transport. As a general 
statistical rule, when the size of the system decreases, 
fluctuations play an increasingly important role. But un- 
der which conditions do they play a role ? and which role 
could they play in nanofluidic transport ? 

There are quite a few studies on this question in rela- 
tion to nanofluidic transport. One may cite the experi- 
mental study of noise in solid state nanopores by Dekker 
et al. |127j . The question of thermal fluctuations was also 
discussed in the context of numerical molecular dynamics 
simulations of osmotic flow through carbon nanotubes by 
Hummer et al. [1281 1129j . Simulations show that - in the 
single flle regime - the flow rate is essentially governed by 
thermal fluctuations rather than hydrodynamics. Hum- 
mer et al. proposed a ID random walk to describe the 
nanofluidic water flow in this stochastic regime. Now, it 
is fair to remark that fluctuations are indeed expected to 
play a key role for single flle transport, for which dynam- 
ics are highly correlated. 

Aside from this situation, various indications of the 
role of thermal fluctuations have been reported, espe- 
cially in studies involving capillary dynamics. One may 
cite the noise effects in the breakup of fluid nanojets by 
Moseler and Landman |130j . the noise assisted spread- 



ing of drops |131j . and the influence of thermal noise in 
thin fllm dewetting |132j in relation to experiments |133j . 
In these different cases, thermal noise modifies quite 
strongly the dynamics, even at relatively large length 
scales (up to typically lOOnm in Ref.[132j. These dif- 
ferent situations can be described within the framework 
of fluctuating hydrodynamics fl34', in which a random 
stress tensor is added to the Navier-Stokes terms. While 
this approach is indeed fruitful in describing the role of 
thermal fluctuations, a general criterion to quantify the 
importance of noise is however lacking. 

Here, to illustrate the potential importance of fluctua- 
tions, we consider a simple situation in which an osmotic 
flow is driven across a single nanopore. The pore has di- 
ameter d and length L and is impermeable to the solute. 
The pressure drop across the membrane is Ap — ksTAc, 
with c the solute concentration. Under the pressure drop 
Ap, a fluid flow is induced, with mean velocity u. As- 



Ap 



The flow 



suming NS equations to hold gives u 
is in the direction of higher solute concentrations. 

Now one may ask the question: could fluctuations yield 
a reverse flow, i. e. against the pressure drop, at least for a 
short time ? This would be clearly a (punctual) violation 
of second law. And what is the minimal size at which 
this may occur ? A lead to such questions is provided 
by the so-called Fluctuation Theorem, which quantifles 
the probability of second law violation by fluctuations. 
This domain is presently a very active, in particular in 
the context of single molecule spectroscopy |135j . The 
theorem states that the probability P{Qr) to find a value 
Qr of the amount of heat dissipated in a time interval r 
satisfies in a non-equilibrium state |136l 11371 1138] 



P{Qr)/P{-Qr) = eMQr/kBT) 



(44) 



In the nanofluidic situation, this theorem would thus 
quantify the probability of a flow against the osmotic 
gradient. Let us compute the order of magnitude for the 
averaged Qr in the stationary state: Qr ~ F x u x t, 
with F the frictional force acting on the nanopore sur- 
face. One has F = TrdLau, with a^i the stress at the 
wall given by ^ dAp/L. To flx ideas, we compute 
the amount of heat dissipated over a time needed for a 
molecule to pass through the whole pore, i.e. t — L/u. 
Altogether this gives 



Qr ^ ApxV 



(45) 



with V = jd'^L the volume of the pore. According to 
Eq. (44 1, violations of the second principle are more 
likely to occur when Qr is of the order of the thermal 
energy. Putting numbers, if we choose Ac = IM (close to 
physiological conditions), then the condition Qr ~ ksT 
is obtained for nanometric volumes, i.e. for a pore with a 
size (diameter, length) typically in the nanometer range. 
This is indeed typically the order of magnitude for the 
size of biological pores. 

The occurence of such 'counter flow' fluctuations has 
actually been observed in MD simulations of osmotic 
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flows through carbon nanotube membranes by Hummer 
et al. |128j . and a stochastic model was proposed to ac- 
count for the water transport across the tube. The above 
analysis thus fixes the limit where fluctuations starts to 
be predominant over the mean behavior. 

Finally we note that one may perform the same esti- 
mate for a different situation of an electric current in- 
duced by an electrostatic potential drop AV^. Assuming 
a bulk conductance, the dissipated heat is estimated un- 
der quite similar terms as Qt ~ elSV x PsV, with pg the 
salt concentration. Again the condition Qt ^ fc^T is 
achieved for nanometric volumes. 

These results point to the crucial role of fluctuations 
in pores with nanometric size. This is the typical scale 
of biological pore: it is then interesting to point out that 
biological pores are working at the edge of second law 
violation (Qr ^ ksT). 

These questions would deserve further experimental 
and theoretical investigations. 



VI. DISCUSSION AND PERSPECTIVES 

In this review, we hope to have convincingly shown 
that nanoscales do indeed play a key role in fluidics. A 
broad panel of length scales ranging from the molecu- 
lar to the micrometric scales leads to a rich ensemble of 
nanofluidic phenomena. While the domain of validity of 
Navier-Stokes equations was shown to extend down to 
the nanometer scale, speciflc transport phenomena show 
up due the complex couplings which build up between 
flow and ionic transport, electrostatics, surface dynam- 
ics, etc. 

We focused in this review on the behavior of fluids at 
the nanoscale. But in doing so we omitted several impor- 
tant topics in the discussions, however strongly connected 
to nanofluidics questions. 

This concern in particular the fleld of transport of 
(macro) molecules, polymers, polylelectrolytes or biolog- 
ical molecules (DNA, RNA) through nanopores. This 
involves either biological nanopores, like the widely stud- 
ied a-haemolysin, or artiflcal solid-state nanopores made 
by ionic drilling of membranes [139] . Starting from the 
pioneering work by Bezrukov et al. and Kasianowicz et 
al. [50l mi [52] , there has been a thorough exploration of 
the translocation mechanisms of macromolecules in tiny 
pores [1401 1141] , with important implications for the un- 
derstanding of biological translocation process and po- 
tentially for single molecule analysis. This domain is now 
rapidly expanding and we refer to the recent review by 
C. Dekker for further reading |139] . 

Another aspect that we left aside is the question of 
nano- and micro- structured surfaces, in line with the re- 
cently developped super-hydrophobic surfaces exhibiting 
the Lotus effect. Superhydrophobic coatings can lead to 
huge slippage effects, with slip lengths in the microm- 
eter range |1SJ [371 11421 1143] , as well as other dynamic 
phenomena, such as hydro-elastic couplings at the super- 



hydrophobic interface |144] . Such surfaces offer the possi- 
bility to considerably enhance the efficiency of transport 
phenomena in particular in the context of slip enhance- 
ment discussed above. For example, as we pointed out in 
the text, a massive ampliflcation, by a factor up to 10^ 
(!), is predicted for diffusio-osmosis on superhydrophobic 
surfaces, which offers the possibility to device efficient 
salt pumps [112] . The underlying coupling mechanisms 
remain however subtle, and this enhancement was shown 
for example to break down for electro-osmotic transport 
|112l I145j . The implications of nano- and micro- struc- 
turation on fluidic transport is still at its infancy and 
remains to be thoroughly explored. It offers the possi- 
bility to couple fluid dynamics over the scales, from the 
nano- to the micro-scales, and even to larger macroscopic 
scales [146 . 

Furthermore, beyond these questions, it is interesting 
to note that nanofluidics offer the possibility to attack 
old questions with new point of views, and new control 
possibilities. We quoted the amazing flow permeabil- 
ity measured by carbon nanotube membranes. While 
this performance remains to be understood - and possi- 
bly investigated at the individual nanotube level this 
opens very promising application in the fleld of energy, in 
particular for portable energy sources. Such membranes 
indeed offers the possibility of a very efficient conver- 
sion of hydrostatic energy to electrical power at small 
scales |1021 I103j . see also |147j . As we have discussed 
above, low (nanoscale) friction at surfaces considerably 
amplifies electrokinetic processes, such as electro-osmosis 
or streaming currents |98l llOlj . The latter phenomenon 
is able to produce an electrical current from a pressure 
drop iSHl [135] • On this basis it was argued [TOU [TU5] 
that the efficiency of this mechano-electric transduction 
raises considerably for slippery walls : reasonable values 
suggest an increase of the energy conversion efficiency 
from a few percents to 30% for slippery surfaces ! While 
the tremendous effect of low surface friction and slip- 
page on charge transport has been recently confirmed 
experimentally jlOlj . its expected impact on mechano- 
electric energy conversion has not received an experimen- 
tal confirmation. Extrapolating with the results obtained 
with the nanotube membranes, this power conversion ef- 
ficiency raises to close to 100% (due to the expected ex- 
tremely weak dissipation) and this amazing prediction 
would lead to a production of electrical power of several 
kW/m2 for a pressure drop of Ibar [102_: this is an im- 
pressive result, which has not been confirmed up to now. 
While such predictions have to be taken with care, they 
strongly suggest to explore the role of nanofluidics in the 
context of energy conversion. Such nanofluidic energy 
conversion devices have a priori the potential to power 
larger scale systems (for example in cars or portable de- 
vices), due to their high power density and low weight 

An alternative fleld of application for the progress 
made in nanofluidic transport is desalination. Indeed, as 
half of the humanity has no immediate access to potable 
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water, desalination of seawater has emerged over the re- 
cent years as an alternative solution to develop fresh wa- 
ter. Reverse osmosis is one of the techniques used in this 
context, for example in the Ashkelon plant in Israel [150]: 
seawater is pushed through a membrane impermeable to 
the salt. Such a process is expensive (typically around 
0.5$/m'^ for the Ashkelon plant in Israel |T5(J'), due - 
among other factors - to the energy required to this op- 
eration. The possibility of reducing this energy by using 
a considerably more permeable membrane should have 
therefore a direct impact on the cost of produced water. 
This is a potential application of membranes made of car- 
bon nanotubes which, as we discussed OISSIES], exhibit 
a permeability 2 to 3 orders of magnitude higher than 
those with micropores, a fully unexpected and still de- 
bated result. Another question for the desalinated water 
is the mineral composition of the produced water which 
- if not controlled - could have a deep, negative, impact 
on agriculture and health. For example the concentra- 
tion of Boron (B), which is toxic to many crops, is high 
in seawater, and requires to be partially eliminated by de- 



salination post-treatments |150j . Reversely, desalination 
removes some ions which are essential to plant growth, 
such as Mg^+ or S04^. There is therefore a potential 
need for a better selectivity of filtration. 

Achieving selectivity similar to biological channels, 
such as aquaporins [HI HH [H], is still out of reach, but 
nanofluidics has a major role to play in this context in 
order to find the key tuning parameters and architecture. 
Such applications may still be years away, but the basic 
principles underlying their operation will develop in the 
years to come. 
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